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1. INTRODUCTION 


1 . 1 Background 

The Stratospheric Aerosol Measurements II (SAM II) and Stratospheric 
Aerosol and Gas Experiment I (SAGE I) satellite experiments are respectively 
the second and third in a series of satellite-borne solar radiometers 
designed to monitor the concentration of aerosols and trace gases in the 
stratosphere. Various other in situ and remote sensing methods have also 
been developed to measure the concentration and variation of stratospheric 
constituents. Except for airborne systems, these other systems provide data 
at specific latitudes and longitudes only, and satellite measurements have 
made available, for the first time, information on a global scale. 

The objective of the SAM II experiment was to obtain vertical profiles 
of the stratospheric aerosol with a vertical resolution of 1 km over the 
latitudinal range obtainable for solar occultation observations from the 
Nimbus-7 solar orbit, approximately 64°-80°N and 64°-80°S. These profiles 
were to be used to study the longitudinal and temporal variations in aerosol 
concentrations. The SAGE I experiment, in addition to measuring 
stratospheric aerosols, also measured the concentration of ozone and nitrogen 
dioxide. In contrast to SAM II, the latitudinal coverage of SAGE I extended 
from about 72°N to 72°S (some seasonal variation occurred within these 
latitude limits). Taken together, the satellites provide almost global 
coverage for the stratospheric aerosol, enabling phenomena such as volcanic 
injection and seasonal variations to be studied in detail. 
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The ultimate objective of the SAM II and SAGE I experiments lies in the use 
of the data to further our knowledge of the processes taking place in the 
stratosphere* Before this can be done, these data must be validated. In 
this process, the data are compared with preexisting models of stratospheric 
aerosols and gaseous constituents, as well as compared with correlative data 
sets, many of which are obtained from specifically designed correlative 
measurement series. The statistical nature of the data has to be examined 
for internal consistency and where SAM II and SAGE I independently make 
measurements at approximately the same time and location these should be 
intercompared. In this report, we describe the results of part of this 
validation process in which the data have been used to study a variety of 
physical, dynamical and chemical processes in the stratosphere. In the 
course of this study, which has yielded a quantity of new scientific 
information, the data have been carefully screened and compared with relevant 
models and data sets. 

1.2 SAM II and SAGE I Satellite Measurements 

The SAM II and SAGE I satellite sensors both consist of sun photometers 
designed to measure atmospheric extinction in certain wavelength pass bands 
during satellite sunrise and sunset (McCormick et al., 1979). SAM II has a 
single pass band centered at 1.0 pm where gaseous absorption is negligible, 
designed to measure extinction due to atmospheric molecular and aerosol 
scattering and absorption. The molecular component of the extinction is 
removed during data processing providing just that due to aerosols. SAGE I 
has four channels, including two centered at 1.0 and 0.45 pm, designated for 
aerosol measurement. The other two channels are at wavelengths of 0.60 pm 
for the measurement of ozone and 0.385 pm for the measurement of nitrogen 
dioxide , respectively. 
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'** The SAM II instrument is mounted, along with other sensors, on board the 
Nimbus 7 satellite* The orbit of this satellite is such that approximately 
14 sunrise and 14 sunset events are observed in one day. The orbit is a high 
noon, sun-synchronous one with all sunrises occurring in the Antarctic region 
and all sunsets in the Arctic region. The exact latitude of the observations 
varies slowly with season between about 64° and 80°, successive events on the 
same day being separated by about 26° in longitude. Nimbus 7 was launched 
October 24, 1978, and continuous SAM II measurements have been made since 
then and the present date. This report includes analyses carried out on the 
data for the period October 1978 to February 1983. SAGE I was launched 
February 18, 1979, on a dedicated AEM-B (Application Explorer Mission B) 
spacecraft into a non-sun synchronous orbit with an inclination of 55°. 
Successive events on the same day are separated by about 24° in longitude, 
and the point of measurement changes slowly in latitude between about 60°N 
and 60°S (the value of latitude extreme depends upon time of year and can be 
as much as 72°N and 72°S) taking 4-6 weeks to complete a cycle. As with SAM 
II, SAGE I initially made observations at both sunrise and sunset, a total of 
30 events per day. Partial failure of the spacecraft power subsystem reduced 
this to sunset events only, about 5 months after launch. Sunset events 
continued to be recorded until final failure of the satellite power 
subsystem, about 34 months after launch. 

As one of these satellites approach the shadow side of the earth 
(spacecraft sunset), the solar radiation reaching it passes through 
successively lower levels in the atmosphere. During this period, the 
radiometer locks onto the sun in azimuth and scans vertically across the 
solar disk. In the course of a sunrise or sunset event, measurements 
covering the whole stratosphere will be made within a period of 20-30 second 
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during which time the sun will be scanned approximately 10 times. The 
atmospheric attenuation is measured along an optical path approximately 300 
km in length through the atmosphere with an instantaneous 0.5 km field of 
view. The data are then inverted to obtain the equivalent vertical 
extinction at each of the radiometer wavelengths (Chu and McCormick, 1979). 
In this inversion procedure, the assumption of horizontal homogeneity of the 
atmospheric constituents is made. In the case of SAGE I, it is necessary to 
separate the contributions to the slant path extinction from each species, as 
well as that due to the molecular atmosphere at each of the four wavelength 
bands. Following this step, the individual vertical profiles for each 
species may be obtained. Separation of contributions from each species 
increases the uncertainties in the data product. An estimate of this 
uncertainty is made along with the inversion and forms part of the final data 
set. 

A significant part of the total post-launch data validation process is 
the operation of field experiments in which correlative data are collected on 
the aerosols and gaseous constituents. These measurements have been made 
from both ground based and high altitude platforms. The former includes 
lidar measurements of the stratosphere aerosol; the latter includes balloon 
and aircraft, in situ particle counters and lidar, rocket and balloon-borne 
ozone measurements and balloon-borne nitrogen dioxide measurements. In the 
present study, the emphasis has been less on the study of individual 
correlative measurements than of the overall data set characteristics. A 
major portion of the entire SAGE I and SAM II data sets have been scanned and 
used to create user-oriented data products. These have included visual 
displays, reference publications and data listings. In addition, the data 
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has been used for global analyses, modeling and related theoretical studies, 
as well as comparison with other long-term data sets. All of these processes 
have provided information on the internal data quality, individual errors and 
relationship to other data sets, as well as unique new information on 
stratospheric behavior. 

1.3 Tasks and Organization of Report 

The Statement of Work, for this contract (NAS1-17032) , specifies eight 
task areas. Work has been carried out under each of these areas and this 
work is reported in sequence in Sections 2 through 7 of this report. The 

title given to each main section is the same as the work statement for Task 
Areas 1 through 6. Task Areas 7 and 8 of the Statement of Work significantly 
overlap other task areas and the description of the work done under these 

two areas has been incorporated with other sections of this report , as 

appropriate. Specifically, the work carried out under Task Areas 7 and 8 is 
reported in the following sections. 

Task Area 7. - Use SAM II and where possible, SAGE satellite data to 

investigate the dynamics of the polar winter stratosphere. 

Work described in Sections 2.8 and 3.1. 

Task Area 8. - Investigate methods of automatically processing and displaying 
SAGE satellite data, such as aerosol optical depth, or suitable map 

projections showing the global variations of the quantity under study. 

Work described in Sections 2.2 and 4.1. 

Specific mention is also made within the body of the report whenever 
work carried out under Task Areas 7 and 8 is described. It should also be 
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noted that although work has been carried out under each task area, the 
division is not necessarily uniform. This division represents the results of 
discussions between the contractors and the NASA Technical Monitor as to the 
work areas requiring major attention. Subsections under each section 
correspond to separate studies within the task area. 

A significant part of the work carried out has been presented at 
scientific meetings and has appeared in scientific journals. Where material 
has been published, the corresponding abstract has been included in the 
Appendix to this report and only a brief summary, together with the 
appropriate references, included in the text. Unpublished material, where it 
is believed to be useful and relevant, has been reproduced in detail in the 
text of this report. 
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2. TASK 1 — INVESTIGATE THE GLOBAL CLIMATOLOGY OF TROPOSPHERIC AND 

STRATOSPHERIC AEROSOLS, AND STRATOSPHERIC 0 3 AND N0 2 INCLUDING SPATIAL 
AND TEMPORAL VARIABILITY 

2.1 Production of Zonal Mean Plots and Tables for SAGE I Data (SAGE I, 
Reference Publication) 

As indicated in Section 1.2, there are 34 months of aerosol, ozone and 
nitrogen dioxide data produced from observations made by the SAGE I satellite 
instrument. This entire data set has been stored, by month, on 34 magnetic 
tapes and archived at the National Space Sciences Data Center, NASA Goddard 
Space Flight Center, Greenbelt, Maryland 20771, where it is available to the 
scientific community. In addition, in order to provide the same community 
with the data in a more easily accessed form, reference publications have 
been developed containing tables and plots of the most significant 
characteristics of the data set. To achieve this goal, a ready-to-use format 
was developed by a small team under the leadership of Dr. M. P. McCormick of 
NASA/LaRC, with assistance from IFAORS and SASC Technologies Inc. Two volumes 
of the product for the SAGE I aerosol data have been produced as NASA 

reference publications (McCormick, 1985, 1986). The publication of the third 
year of aerosol results is expected to follow soon. 

2.2 Production of SAGE I Aerosol Optical Depth Maps and Zonal Mean Extinction 

Plots (includes work carried out under Task Area 8) 

One of the more useful quantitative measures of the stratospheric 
aerosol concentration that may be derived from the SAM II and SAGE I data 
sets is the optical depth at a wavelength of 1 pm. In order to study this 
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quantity on a global scale, maps of stratospheric optical depth have been 
prepared for the complete SAGE I data set. The original production of these 
maps was manual and both color and black and white versions were prepared for 
some periods of data. Some of these maps have been used in publications and 
an example (taken from Kent and McCormick, 1984), is reproduced in Figure 
2.1(a). A comparison data product, showing the zonal mean extinction as a 
function of altitude and latitude, has also been prepared and is shown in 
Figure 2. 1(b) . 

Manual production of these types of display is not only lengthy and 
tedious, but can introduce subjective error into the interpolation of the 
data points. For these reasons, software has been written to produce the 
optical depth maps automatically by computer. This problem is not simple, as 
the SAGE data points are not distributed uniformly across the globe, but over 
a five or six-week period will scan the earth f s surface in latitude from one 
extreme to another. An example of this work is shown in Fig. 2.2. In 
addition, data points may be missing, necessitating an interpolation 
procedure. The software to plot the required map has been divided into three 
parts as described below. 

(a) Data Selection and Interpolation 

Independent maps are prepared for sunrise and sunset events (only sunset 
after June 1979). The maps show the optical depth for a sunrise or sunset 
latitude "sweep,” such as that in Figure 2.2. Sunrise or sunset data is 
selected that lies between two chosen dates that define the ends of the 
sweep. Individual events are sequenced from one end of the sweep to the 
other, the latitudes being examined to ensure that the event sequence moves 
steadily in one latitudinal direction. Following this, interpolation of data 


8 


ORIG’NAL PAGE TS 

OF POOR QUALITY 


m: 




9^+1 


SHADING ° d P e T ^ L SHADING ° d£ ^ L 


cm < 1 . 0*10 , 
r — i 10 - 20 * 10 '^ 
C— 1 2.0-3 Ox io - 5 


P:':^ya 30-40*10 

■ 40-5 0*10' 

mm > 5 o*io“ 3 


90 1 I I I L 1 . J 1 I I I L.-.l 

-IBO-ITOieO^SO-MO-tM-GO-ltO-IOO-gO-SO^O-eO-SO^O -30-20 -10 0 K5 20 30 4050&070 BO 90 OO 1C !20 130 *0 *0 


LONGITUDE 


Figure 2.1(a) 


SAGE map of the stratospheric optical depth, July 2 to 
August 10, 1981. Note the concentration of volcanic material 
to the north, due to the eruption of Alaid and along the 
equator, due to the eruptions of Pagan and Ulawun. 
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Figure 2.1(b). SAGE measurements of zonal mean extinction July 2 to August 10 
1981. The baseline below the shading is at an altitude 2 km 
above the tropopause. Dashed vertical lines show the mean 
latitude for each day's measurements. Note the concentrations 
of volcanic material due to the eruptions of Ulawun (5°S), 
Pagan (18°N), and Alaid (51°N). 


gaps is carried out. Two classes of data loss may be identified. The first 
consists of those unavoidable occasions in which the satellite is in full 
sunlight and no data is obtainable. No interpolation is attempted on these 
occasions. The second class consists of those occasions on which system 
malfunctions (transmission drop-out, noise, etc.) have resulted in the loss 
of one or more data events. A computer routine identifies these events and, 
provided the loss does not involve too many geographically adjacent 
measurement locations, it carries out an interpolation between the 
neighboring points in a suitable manner so as to fill in the missing data. 
In this way, a typical latitude sweep (e.g., August 28-September 25, 1980) 

consisting of about four hundred observations (between latitudes of 66°S and 
60°N) with about thirty missing observations is converted into a complete 
data sequence with no gaps. 

(b) Interpolation to a Rectangular Grid 

Even given a complete data sequence, problems still exist in its 
mapping. As noted above, the data points are non-rectangularly distributed, 
the latitude and longitudinal separations are different and the latitude 
separation is highly dependent on location on the earth's surface. Thus, at 
the equator, points at approximately the same longitude on successive days 
will be separated by about 5° in latitude. At high latitudes, this 
separation may be less than 1°. A decision was made to interpolate all data 
to a 1° latitude, longitude grid, this being a suitable base for plotting 
onto a mercator projection of the earth's surface. A computer routine 
divides the data sequence into skew quadrilateral areas on the earth's 
surface. The optical depth values on a rectangular 1° latitude-longitude 
grid within each quadrilaterial are interpolated by a three-dimensional 
projection of the quadrilaterial onto a rectangle, followed by a standard 
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bivariate interpolation routine. The resultant data matrix, consisting of 
values of optical depth at 1° latitude-longitude intervals is output and 
passed to the plotting routine. 

(c) Production of global map 

The plotting routine produces three maps. 

(1) The first shows the measurement locations, an example has already 
been given in Figure 2.2. 

(2) The second map is similar to the first with the optical depth 
values superimposed onto the measurement locations. An example is 
given in Figure 2.3. 

The third is a shaded contour map prepared from the one degree 
interpolated optical depth grid described previously. Six shading 

levels are used. Considerable effort has been employed in creating a 
suitable shading scale with the required resolution. All maps are 
drawn 9" x 18" making 1° of latitude or longitude equal to 0.05". This 
represents the limits of resolution of the figure; when the figure is 
reduced to page size or smaller, the associated discontinuities become 
extremely small. The boundaries of the shaded regions are enhanced by 
the addition of contour lines. An example of this plot, showing the 
same data as in Figures 2.2. and 2.3, is reproduced in Figure 2.4. 

Plotting time, particularly for Figure 2.3, is lengthy and can be 
expensive if much of the map is covered by high density shading. Plots of a 
quality suitable for scientific use may be made on the Varian plotter at the 
NASA-Langley Computing Center for which the charges are minimal. Plots of 
publication quality must be produced on the Calcomp plotter at considerably 
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Figure 2.2. Positions of SAGE observations between May 11, 1980 and June 23, 1980. 

Positions are shown by plus signs. Missing data points, for which 
interpolation has been carried out, are shown by circled plus signs. 
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Figure 2.3. Optical depth values (from 2 km above the tropopause) for the observation 
points shown in Figure 2.2. 






more cost. It may be noted that the output matrix file from stage (2) above may 
be used as input to other plotting devices to produce alternative color 
plots, vugraphs or black and white plots at less cost. 

Routine production of the SAGE I optical depth maps has been carried out 
for all useful data between launch in February 1979 and December 1980. During 
this routine production, the computer codes were refined, particularly in 
terms of handling data gaps. These occurred in the SAGE I data sequence 
because of (a) sunlit periods, (b) telemetry problems, and (c) spacecraft 
power problems. Gaps due to sunlit periods are large and present no major 
plotting problems. Gaps due to telemetry and spacecraft power supply 
problems are erratic and of variable size. As noted earlier, the code 

contains an interpolation routine which will interpolate across data gaps 
that are below a certain size. This routine is capable of handling most of 
the data except during the period May-October 1979, when the data is very 
irregular due to partial failure in the spacecraft power supply. For this 
and other similar periods, the code was modified to plot only complete data 
segments after normal interpolation has been attempted. Up to ten separate 
data segments per sweep may be plotted although this has never been necessary 
in practice. 

In addition to data gaps, problems have been encountered due to poor 
satellite ephemeris coordinates. Occasionally, the latitude (or longitude) 
of an observation is out of sequence. The code will detect and correct 
single errors of this type. If, however, as occasionally occurs, such an 
error is adjacent to a data gap or two such errors occur in sequence, then it 
is necessary to correct the latitude and longitude on the original data 
file. A list of optical depth maps prepared is shown in Table 2.1 which 
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TABLE 2 . 1 . Stratospheric Aerosol Optical Depth Maps 
(February 21 , 1979 - April 8 , 1980 ). 
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TABLE 2.1. Stratospheric Aerosol Optical Depth Haps (cont'd) 
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indicates the dates of each sweep , the event type , sweep direction and the 
plotting device used (these maps have been delivered to NASA). 

2.3 Ozone Climatology 

The entire three years' (34 months) SAGE ozone measurements have been 
used to derive the seasonal meridional ozone mixing ratio distributions. The 
results are tabulated in Tables 2.2 to 2.5 and graphically displayed in 
Figures 2.5a to 2.5d. In general, the ozone mixing ratio shows a maximum 
above the equator at S&32 km. The center of this maximum appears to follow 
the seasonal shift of the sun. We may also note that, in the winter high 
latitude lower stratosphere (between 25 and 35 km), the ozone shows a lower 
value than that in the corresponding region in the summer hemisphere. This 
occurrence could be the result of large-scale wave disturbances (planetary 
waves) in high latitudes during the winter. Planetary waves are believed to 
be able to transport trace species and heat to the polar region from lower 
latitudes. 

The corresponding standard deviations (expressed in percentage of the 
mean) of the three-year seasonal means are given in Figures 2.6a to 2.6d. 
Generally, large values (>20%) of the standard deviation exist in regions 
above 45 km and below 20 km, especially in the low latitudes (<30°) in the 
lower stratosphere. Perhaps, these large values are partly due to large 
signal/noise ratios for the instrument in these altitude ranges. In the 
regions between 20 km and 45 km, the values are generally less than 20%. In 
fact, values less than 10% are found in the stratosphere between 23 and 38 
km, except in the high latitudes. Since in this region the uncertainty of 
the retrieved 0^ profiles is less than 10%, the larger values (>10%) in the 
high latitudes are indicative of the ozone variability in this region. Of 
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TABLE 2.4. Zonally Averaged, Meridional Ozone Mixing Ratio (PPMV) for Northern Fall-Southern 
Spring. Derived from the Three-Year SAGE Ozone Measurements. 
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TABLE 2.5. Zonally Averaged, Meridional Ozone Mixing Ratio (PPMV) for Northern Winter Southern 
Summer . Derived from the Three-Year SAGE Ozone Measurements. 
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Figure 2.6. Seasonal meridional ozone mixing ratio — percentage standard deviation from mean. Derived 
from 3-year SAGE ozone measurements (a) March to May, (b) June to August, (c) September 
to November, and (d) December to February. 





particular interest is the extinsionof the the 10% contour line from winter high 
latitudes to about +10° centered at altitude 32 km (Figures 2.6b and 2.6d). There are 
reasons to believe that this feature is primarily a result of planetary wave 
activity during the winter. Since these wave activities are prevalent 

through about the end of March (Labitzke, 1982, J. Meteo. Soc^, Japan , 

124-139) relatively high values of standard deviation are also expected to be 
shown in high latitudes during the period of spring equinox. This feature is 
indeed noticeable in Figures 2.6a and 2.6c. 

For the reason that the ozone column plays an important role in uv solar 
radiation attenuation, extensive investigations have beem made of the 
variation of total ozone on a global scale. To examine the seasonal 

variation of the zonally averaged total ozone distribution, three-year ozone 

* 

profiles derived from the SAGE instrument have been used for this particular 
purpose. In this analysis, the ozone column is integrated from 2 km above 
the tropopause up to 55 km. The 34-month SAGE ozone profiles are grouped 
seasonally. The ozone column is then determiend for each of the profiles. A 
zonal average is then calculated for the 15 10° latitude bins from 75°S to 
75°N. The associated standard deviations of the means have also been 
derived. These deviations may indicate the variabilities of the ozone column 
at different latitudes during different seasons. The four seasonal means are 
shown in Figure 2.7. In general, large seasonal variations of ozone columns 
exist in regions outside the tropics. In both hemispheres, Figure 2.7 shows 
that the mean ozone column reaches its lowest value in the Fall. A rapid 
increase occurs from the Fall to Winter at high latitudes. This increase 
seems to extend to Spring in the Northern Hemisphere. In the Southern 
Hemisphere, the variations are similar to those in the Northern Hemisphere, 
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except poleward of 50°S. At the latter latitudes, the ozone column decreases, 
rather than increases from Winter to Spring. 

The standard deviations of the seasonally averaged zonal mean ozone 
column are given in Figures 2.8a to 2.8d. The remarkable features are the 
larger values of the deviation during the Winter and Spring, and the 
relatively low values during the Fall in high latitudes. In general, the 
values for the deviation in the Northern Hemisphere are generally greater 
than those in the Southern Hemisphere during the same seasons. These large 
deviations at the high latitudes during Winter and Spring, especially in the 
Northern Hemisphere, are believed to be linked with the high latitude dynamic 
activity during these particular seasons. It is well known that the 
stratosphere is characterized by large-scale disturbances (planetary waves) 
as the circulation pattern changes from Winter to Summer. At low latitudes, 
the standard deviation is much smaller compared to those at high latitudes. 
Nevertheless, Figure 2.8 shows a mild local maximum of the deviation at the 
tropics and two local minima near +20 degrees latitude. 

2.4 Identification of Calibration Orbits 

The SAM II satellite system is mounted, along with other sensors, aboard 
the Nimbus-7 spacecraft. Periodic earth radiation budget instrument 
calibrations, carried out on board Nimbus-7, result in perturbations to the 
spacecraft attitude control/pointing which interfere with the data being 
measured by SAM II. These calibrations are carried out at approximately 
24-day intervals and perturb only a small fraction of the SAM II measurement 
orbits. Despite this, the interference, which simulates a high extinction at 
altitudes above the aerosol layer maximum, is sufficiently large to seriously 
affect average values that include the perturbed data. 
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In order to identify SAM II events that are affected by calibration, a 
computer search was made for the apparent high extinction that appears on the 
SAM II extinction profile. The computer program integrates the SAM II 
extinction profiles over the altitude range 26-30 km and has been applied to 
all data obtained between 1978 and 1981. For later years, which contained 
high altitude extinction enhancements as a result of the eruption of El 
Chichon, the integration range was shifted to 31-35 km. A list of orbit 
numbers for those orbits then judged to be affected by calibration are listed 
in Table 2.6. In this table, asterisks . mark those orbits for which no clear 
conclusion was possible. Data that were not available on the microfiche, but 
which showed an obviously high extinction, are shown in Table 2.7. 

2.5 Production of a 16 mm Movie Film of SAM II Data 

In order to observe the temporal variation of the aerosol extinction 
coefficient, especially the generation, movement and disappearance of polar 
stratospheric clouds, a computer program was written to generate a 16 mm 
movie film from the SAM II data set. In this movie, the aerosol extinction 
contours are shown within a height versus longitude frame. The height range 
is fixed from 5 km to 30 km, while the longitude range is fixed from -180° to 
180°. SAM II sunrise observations occur only in the Southern Hemisphere while 
sunset observations occur in the Northern Hemisphere. In order to show the 
aerosol behavior in each hemisphere independently, the sunrise and sunset 
events are separated into two different sequences. Each frame of the movie 
for either hemisphere is generated using a single day’s sunrise or sunset 
data. The date and the latitude of the first observation on each day are 
shown at the top of each frame. In addition, the position of the tropopause 
is shown by a solid line with dots on it. 



TABLE 2.6. Calibration Orbits Occurring in SAM II Data. 
November 1978 - February 1983. 


Year 

Orbit Number 

Month /Day 

SAM II Week 

ISUN 

1978 

A 90 

N o v 2 B 

nnA 

VV V 

r\ 

\j 


797 

Dec 20 

009 

i 

1979 

1552 

Feb 13 

017 

0 


1773 

Mar 01 

019 

0 


2105 

Mar 25 

023 

0 


2105 

Mar 25 

023 

1 


2568 

Apr 28 

027 

0 


2568 

Apr 28 

027 

1 


2928 

May 24 

031 

0 


3259 

Jun 16 

034 

0 


3259 

Jun 16 

034 

1 


3591 

Jul 11 

038 

0 


3591 

Jul 11 

038 

1 


3923 

Aug 04 

041 

0 


4254 

Aug 28 

045 

0 


4586 

Sep 12 

048 

0 


4586 

Sep 12 

048 

1 


5097 

Oct 28 

054 

1 


5249 

Nov 08 

055 

1 


5560* 

Nov 30 

058 

1 


5581 

Dec 01 

059 

0 


5915 

Dec 26 

062 

0 

1980 

6244 

Jan 19 

065 

0 


6576 

Feb 12 

069 

0 


7239 

Mar 31 

076 

0 


7571 

Apr 24 

079 

0 


7903 

May 18 

083 

0 


8567 

Jul 05 

089 

0 
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TABLE 2.6 (Cont'd) 


Year Orbit Number Month /Day SAM II Week I SUN 


1980 (cont'd) 8567 Jul 

8897 Jul 

9561 Sep 

9892 Oct 

10225 Nov 

10887 Dec 

1981 11219 Jan 

11549 Feb 

11549 Feb 

11882 Mar 

11882 Mar 

12214 Mar 

12878 May 

13210 Jun 

13873 Jul 

14205 Aug 

14536 Sep 

14868 Oct 

15533 Nov 

15863 Dec 

1982 16195 Jan 

16195 Jan 

16858 Feb 

16858 Feb 

17190 Mar 

17190 Mar 

17521* Apr 

17853 May 

18185 Jun 

18185 Jun 

18517 Jun 

18517 Jun 


05 

089 

1 

29 

093 

1 

15 

100 

0 

09 

103 

0 

02 

107 

1 

02 

113 

1 

13 

117 

0 

05 

120 

0 

05 

120 

1 

02 

124 

0 

02 

124 

1 

26 

127 

0 

13 

134 

0 

06 

137 

1 

24 

144 

1 

17 

148 

0 

01 

151 

0 

04 

155 

0 

21 

161 

0 

15 

165 

0 

08 

168 

0 

08 

168 

1 

25 

175 

0 

25 

175 

1 

21 

179 

0 

21 

179 

1 

14 

182 

1 

08 

185 

1 

10 

189 

0 

10 

189 

1 

25 

192 

0 

25 

192 

1 
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TABLE 2.6 (Cont'd) 


Year 

Orbit Number 

Month /Day 

SAM II Week 

ISUN 

1982 (cont'd) 

18848 

Jul 19 

196 

0 


18848 

Jul 19 

196 

1 


19512 

Sep 05 

203 

0 


19512 

Sep 05 

203 

1 


19844 

Sep 29 

206 

0 


i no/./. 

J. 

Oa« on 

206 

1 

X 

1983 

21834 

Feb 20 

227 

1 


TABLE 2.7. Orbits Showing High Extinction Above 26 km for which 
No Microfiche is Available 


Year 

Orbit Number 

Month /Day 

SAM II Week 

ISUN 

1979 

5130 

Oct 30 

054 

1 


5133 

Oct 30 

054 

1 

1980 

6033 

Jan 03 

063 

1 


6034 

Jan 03 

063 

1 

1982 

21148 

Dec 31 

219 

1 
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Two movies were generated from the developed software. One movie shows 
polar stratospheric clouds in the northern winter of 1978 (sunset 
observations only). The other movie shows polar stratospheric clouds in the 
Southern Hemisphere winter of 1979 (sunrise observations only). The 
development of polar stratospheric clouds and their apparent vertical 
movement associated with the vertical variation of the tropopause height are 
clearly shown in these films (these films have been delivered to NASA) . 

A contour plot of SAM 11 extinction data used in the generation of a 
16mm film is shown in Figure 2.9. 

2.6 Comparison of SAM II and SAGE I Data 

The SAM II and SAGE I satellite experiments are each designed to measure 
stratospheric aerosol extinction coefficients using the technique of solar 
occultation or limb extinction. Both SAM II and SAGE I have an aerosol 
measurement channel centered at 1.0 pm wavelength, and have a similar optical 
bandpass and instrument field of view. There are occasions when the 
measurement locations of the two satellites are nearly coincident, providing 
opportunities for a measurement comparison. 

In this study, the aerosol extinction profile and daily contour plots 
for the coincident events in 1979 have been compared. The first approach was 
to compare individual SAM II and SAGE measurement profiles obtained on the 
same day and within 1° in both longitude and latitude. It was found that 
there are thirteen profiles in 1979 that satisfy this criterion. The second 
approach was to compare the daily aerosol extinction contour plots for the 
two satellites. Again, comparisons were made when the difference between the 
latitude coverages of SAM II and SAGE I was 1° or less. It was found that 
there were eighteen daily contour plots in 1979 that satisfy this criterion. 
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Figure 2.9. Contour plot of SAM II extinction data used in the generation of 16mm film. 






After carefully comparing the individual profiles of aerosol extinction and 
daily contour plots obtained by SAM II and SAGE I for occasions when the 
measurements were nearly coincident in time and space, we concluded that both 

profiles and daily contour plots agreed within the uncertainties associated 
with each experiment. Thus, because SAM II profiles are in good agreement 
with results obtained using more conventional remote and in situ sensors, the 
results of this study validate the data obtained with the SAGE I, 1.0 pm 
sensor channel. 

This work has been published in the Journal of Geophysical Research (Yue 
et al., 1984) and the Abstract to this publication is reproduced as Appendix 
A. 

2.7 Polar Stratospheric Cloud Studies 

Polar stratospheric clouds (PSC's) appear on the SAM II data as 
enhancements in the 1 pm aerosol extinction at altitudes between zero and ten 
kilometers above the tropopause; a few sightings of these clouds have also 
been found in the SAGE I data set. PSC's occur only at high latitude, in 
winter, when the ambient stratospheric temperature falls below about 195K. 
Their characteristics, as seen on the satellite data, and theories of their 
formation have formed the basis of several publications by members of the 
Aerosol Research Branch, LaRC (McCormick et al., 1981, 1982, 1985; Steele et 

al., 1983, Hamill and McMaster, 1984). As part of the work carried out under 
this contract, a study was made of the occurrence statistics for PSC's and 
this work has formed the subject of a NASA Contractor Report (Farrukh 1985). 
Data for both hemispheres were studied for the period 1978-1982. The 
following features were examined in particular: 
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1. The increase in extinction (for a PSC) relative to the extinction of 
the background (non PSC) aerosol. 

2. The altitude variation of minimum stratospheric temperatures. 

3. The longitudinal distribution of PSCs. 

4. The longitudinal distribution of minimum stratospheric temperatures. 

5. The probability of occurrence of PSCs as a function of minimum 
stratospheric temperature. 

Principal findings from this analysis are summarized in the abstract to 
the report cited above (Farrukh 1985) which is reproduced as Appendix B. Part 
of the material has also been included in a further publication (McCormick et 
al, 1985). The abstract for the publication is included as Appendix C. 

2.8 Study of Aerosol Concentrations Near the North Polar Vortex (includes work 
carried out under Task Area 7) 

Airborne lidar measurements made in the arctic region in January 1983 
showed strong gradients in aerosol concentration to exist across the region 
of the polar night jet stream. Similar differences were observed in the SAM 
II data taken at the same time as the lidar flights, indicating that for 
altitudes between 18 and 26 km, the aerosol concentrations within the polar 
vortex were about one order of magnitude lower than that outside. This 

analysis has been published (McCormick et al., 1983) and the abstract is 

reproduced as Appendix D. The observations just outlined were made shortly 
after the eruption of El Chichon had injected a very large quantity of fresh 

aerosol into the stratosphere. It was not clear whether the aerosol 

concentration variations observed were entirely due to the action of the 
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north polar vortex or if there may have been some effect of inhomogeneities 
in concentration remaining from the volcanic eruption. Partly because of 
this uncertainty and also because of a need to investigate the underlying 
mechanisms, an indepth study was made of the aerosol extinction in high 
latitudes in the northern winter of 1978-1979, as measured by SAM II. These 
studies confirmed that the concentration changes were produced as the result 
of dynamical changes associated with the north polar vortex. Subsidence was 
observed to occur within the vortex, bringing air down from higher altitudes 
containing a lower aerosol concentration and mixing out air with a high 
aerosol concentration through the base of the vortex. Rates of subsidence 
were established and a base for the vortex was determined at an altitude of 
about 14 km. This work has been published in the Journal of the Atmospheric 
Sciences (Kent et al., 1985a) and the abstract is reproduced as Appendix E. 

2.9 Tropospheric Aerosol Climatology 

As part of a program of work carried out under NASA Contract NAS8-35594, 
a study was made of aerosol extinction values measured by SAM II and SAGE I 
for altitudes below the tropopause (Kent et al., 1985b). The SAM II and SAGE 
I satellites were designed for stratospheric use but it was found that, in 
the absence of high altitude cloud, tropospheric measurements are not only 
possible, but are available from a significant fraction of the total number 
of measurement opportunities. In Che previous work, SAM II data for October 
1978 - October 1979 and SAGE I, 1 pm data for October 1979 - November 1981 
were studied in terms of the probability distribution of extinction values at 
different altitudes and latitudes. The 502 probability level or median 
extinction was chosen as a measure of central tendency for the extinction 
values. The main results to emerge from this analysis were; 


1. Median extinction values in the Northern Hemisphere for altitudes 
between 5 and 10 km were one-half to one order of magnitude greater 
than in the Southern Hemisphere. 

2. A seasonal variation by a factor of 1.5 - 2 was observed in both 
hemispheres with maximum extinction being observed in local Spring and 
Summer • 

3. Following major volcanic eruptions, a long-lived enhancement of the 
aerosol extinction was observed for altitudes above 5 km. 

Further analysis of the tropospheric data from SAGE II, as well as SAM 
II and SAGE I, has continued under Contract NAS 1-18252 and also under the 
present contract. The work under the current contract has concentrated upon 

the SAM II data, of which only one year (October 1978-October 1979) was 

previously studied. The emphasis of the newer work has been to use the 
aerosol/molecular extinction ratio as the parameter to be studied, rather 
than the extinction itself. The former is a quantity which, to a first order 
of approximation, is conserved under transport and thus more suitable for use 
in a region where both horizontal transport and vertical convection can 
occur. In addition, instead of using the median as a description of aerosol 
properties, the probability distribution has been examined for modes (most 
probable values) in the extinction ratio. The variation of these with 

altitude and season has then been studied. In both present and past 

analyses, we have tended to use a logarithmic scale to describe the aerosol 
extinction or extinction ratio. This choice has been necessitated by the 
large dynamic ranges for these variables found in the tropospheric data. 

Figure 2.10. shows a typical extipctipn ratio probability distribution 
for SAM II data. This histogram shows the frequency of occurrence of various 

3 ® 




AEROSOL/MOLECULAR EXTINCTION RATIO 


Figure 2,10. Histogram showing frequency of occurrence of different 
values of aerosol/molecular extinction ratio* SAM II 
data, November 19 - December 16, 1979, Northern hemisphere. 
Data are for an altitude of 7 km. 
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values of extinction ratio for one month of SAM II data taken in the Northern 


Hemisphere. The data is for an altitude of 7 km and was taken between 
November 19 and December 16, 1979. The extinction ratio is shown on a 
logarithmic scale, each interval in the histogram corresponding to a change 
of 0.1 in the logarithm of the extinction ratio. The histogram clearly shows 
two maximas, one at an extinction ratio of about 0.6 and the other at a value 
of about 22. The former is believed to be associated with the tropospheric 
aerosol but the latter, which lies at the threshold of measurement of the 
satellite instrument (at this altitude) is believed to be associated with 

subvisible clouds. Support for this interpretation is given by the 

multiwavelength SAGE I and SAGE II data which show similar modes. For data 
from these two satellites, it is possible to calculate the ratio of the 
extinction at a shorter wavelength (0.45 pm for SAGE I or 0.525 pm for SAGE 
II) to that at 1 pm, for the same altitude and the same event. It is found 
that such ratios, for the mode just identified as aerosol, typically have 
values around 2 or 3. For the other mode, this ratio is normally around 
unity. The latter is to be expected for larger cloud particles, lending 

credence to the interpretation just given. 

It is possible to study the manner in which the tropospheric aerosol 

modes, as just described, vary with both altitude and season in either 
hemisphere. The results of this study are shown in Figure 2.11 in the form 
of a shaded contour plot. The time period covered in the plot is October 
1979 - September 1981, the data being divided into four-week time segments. 
The altitude range is from 5 to 15 km, with a vertical resolution of 1 km. 
The horizontal bars on the diagram show the mean altitude of the tropopause 
for each four-week period. The interval 5-15 km was chosen to show behavior 
in the lower stratosphere, as well as that in the upper troposphere. Some 
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tropospheric data below 5 km is available but the quantity is low and 
statistical fluctuations in the mode values are large; for this reason, it 
has not been plotted. It may also be noted that Figure 2.11 shows some data 
missing above 5 km; this occurrence is common after volcanic injection into 
the upper troposphere and lower stratosphere has raised the extinction to 
such a level that the signal received from the sun at the satellite falls 
below the threshold of detection of the instrument. The quantity plotted in 
the figure is the aerosol/molecular extinction ratio at the aerosol mode as 
shown in Figure 2.10. Six levels of extinction ratio are shown covering most 
of the values encountered in the non-volcanic situation. 

Figure 2.11(a) shows data for the Antarctic Region. In the stratosphere, 
the seasonal cycle is clearly evident with extinction maxima occuring in 
June-August in each year (associated with the formation of polar 
stratospheric clouds). There is also a gradual increase in the extinction 

ratio between 1979 and 1981, most likely associated with the eruption of 
Sierra Negra in late 1979. In the troposphere, there is a high degree of 
vertical uniformity in extinction ratio. There is also a weak seasonal cycle 
with minimum values in the April-July time frame. Figure 2.11(b) shows the 
variation of the same quantity for the Northern Hemisphere. In general, 
values are higher than in the Southern Hemisphere even before the eruptions 
of St. Helens in 1980 and Alaid in 1981 produced very high extinction 
ratios. There is also a strong seasonal cycle with maxima in the summer 
period that is partially masked by the volcanic effects in 1980-1981. Summer 
1979 and late Spring 1980, prior to the May eruption of St. Helens, show 
uniformly high values of extinction ratio up to the tropopause indicating a 
high degree of vertical mixing. The winters of 1978-1979, 1979-1980, and 
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1980-1981 show low extinction ratios (as in the Antarctic Region), even in 
1980-1981 when the volcanic effects might be expected to be important. 

In summary 

1. Extinction ratio values are higher in the Northern Hemisphere. 

2. There is a seasonal cycle in both hemispheres with maxima in local 
Spring-Summer and minima in local fall-winter. 

3. There is a high degree of vertical uniformity in the extinction ratio 
for altitudes between 5 km and the tropopause. 

4. Volcanic effects are strong in the upper troposphere, as well as the 
stratosphere but do not totally mask the seasonal cycle. 

2. 10 The Aging Process of Aerosol Particles in the Stratosphere 

In situ measurements of stratospheric aerosols at different latitudes 
have shown that large numbers of small particles exist in the tropical zone, 
and that the aerosol layer extended from higher altitudes near the equator to 
lower ones at higher latitudes. This behavior indicates that the upper 
troposphere in tropical regions may be a source of condensation nuclei in the 
stratosphere and that these particles mature as they move to higher altitudes 
and latitudes. 

In this study, the aerosol extinctions measured at 0.45 pm and 1.0 pm 
from the SAGE I data set were utilized to deduce the aerosol size 
distribution (See Section 7.1. Retrieval of Aerosol Properties from Aerosol 
Extinction Ratios at 0.45 pm and 1.0 pm). In order to find out the 
latitudinal and altitudinal variation of the aerosol size distribution, the 
atmosphere was divided into a grid consisting of 1 km intervals in altitude 
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and 5° intervals in latitude • The altitudes under consideration ranged from 
10 km to 25 km and the latitudes from 55°S to 50°N. The 25 km upper limit for 
altitude was necessary since 0.45 pm extinction values at altitudes above 25 
km have much larger error bars. For each gridpoint, the average value of the 
ratio of extinctions at the two wavelengths for the same month of SAGE I data 
was calculated. This value was then converted to mode radius in the assumed 
lognormal size distribution of aerosol particles. The variation of mode 
radius with latitude and altitude was then studied. Although only March 1979 
and April 1979 data were analyzed, our results support the view that aerosols 
grow as they rise through the stratosphere as predicted by the modeling 
results of Turco et al., (1979) and concluded by different experimentalists, 
including Bigg (1976), Farlow ( 1979), and Oberbeck et al. , (1981). This work 
has been published in Geophysics Research Letters (Yue and Deepak, Geophys . 
Res . Lett . ,11, 999-1002, 1984). The abstract is included as Appendix F. 


44 


3. TASK 2 — TO INVESTIGATE THE CORRELATION OF STRATOSPHERIC AEROSOLS, 0 3 AND 

N0 2 , WITH EACH OTHER AND VARIOUS METEOROLOGICAL PARAMETERS 

3.1 Correlation Between Aerosol Optical Properties and 30 mb Pressure Height 

(includes work carried out under Task Area 7) 

3.1.1 Introduction 

Aerosol optical depths at 1 jam wavelength above 18 km, derived from the 
Stratospheric Aerosol Measurement (SAM II) on February I, 1983, were shown to 
be lower by approximately one order of magnitude within the polar vortex than 
those outside (McCormick et al., 1983). This result indicates the effect of 
winter circulation on the distribution of aerosols during the winter season. 
The downward motion^ in the polar region during the development stage of the 
circumpolar vortex, can be the major reason for this feature since the 
aerosol extinction decreases with altitude. Furthermore, the relatively low 
aerosol optical depth is also evident during periods of stratospheric 
warmings when the circumpolar vortex is substantially disturbed by amplified 
planetary waves, (Wang and McCormick, 1985; see also Section 5.2). The reason 
for this is that the air mass within the low pressure system(s) is largely 
able to preserve its polar characteristics as the planetary waves intensify. 
Cross-contour mixing is rather inefficient due to the prevailing geostophic 
balance, (Wang and McCormick, 1985; See Section 5.2). Aerosol extinction data 
from both SAM II and SAGE for four selected days in different years taken in 
the high latitude Northern Hemisphere, were used in that analysis. (The 
results obtained have been further substantiated using SAGE 0^ and N0 2 
data.) The results of Wang and McCormick (1985), suggest that a significant 


1. The development of the downward motion and the polar vortex has been 
discussed in the model analyses by Leovy (1964) and Holton and Wehrbein 
(1980), and also in the texts by Craig (1965) and Holton (1975). 
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correlation exists between the aerosol optical depth (above 20 km) and the 
height at 30 mb pressure level. The work presented here is devoted to the 
investigation of this correlation. In this study, we have used almost all 
the SAM II measurements during periods when the stratospheric planetary waves 
are most active. 

3.1.2 Approach 

The Arctic extinction profiles from January to March of the first four 
years of SAM II measurements ( 1979 — 1982) have been used to study the linear 
correlation between aerosol optical depth and 30 mb height. The extinction 
profiles for the same period in 1983 are found showing enhanced aerosol 
loading due to the effects of the El Chichon volcanic eruption. For this 
reason, the SAM II measurements made in 1983 have not been included in this 
study. In addition, large variations in extinction associated with polar 
stratospheric clouds have been reported (McCormick et al., 1982). For our 
purpose, it was necessary to filter out those special events by including 
only profiles above 14 km for which the SAM II extinction ratio is less than 
4. We have also investigated the results obtained when using extinction 
ratios of 3.5 and 4.5 as the criterion. The change in the calculated results 
was found to be less than 2% using these three alternatives. 

3.1.3 Results and discussion 

The analysis of the linear correlation between aerosol optical depth 
(above 20 km) and the height of the 30 mb pressure level for January-March 
1979 is displayed in Figure 3.1. In the derivation of the correlation 
coefficient, the 30 mb height has been taken as the dependent variable. 
There are 1036 data pairs in this 1979 case. The constant (a) and slope (b) 

3 

for a linear least-square fit were found to be 22.1 and 2.69 x 10 , 
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1979 

JANUARY-MARCH 

LINEAR CORRELATION, y = o + bx 
a = 22.l, b=2.69xl0 3 
LINEAR CORRELATION C0EFF = .834 



REROSOl OPTICAL DEPTH. 1 .E-3 


Figure 3.1. Linear correlation between SAM II aerosol optical depth (above 
20 km) and associated 30 mb height for data from January 
to March 1979. 
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respectively. The computed linear correlation coefficient is 0.834. This 
result indicates a significant positive correlation between the aerosol 
optical depth and the 30 mb height. This relationship is also illustrated in 
Figure 3.2 which shows 3-D projective plots of the aerosol optical depth 
(Figure 3.2a) and the 30 mb height (Figure 3.2b) as functions of longitude 
and time (day). The period for these plots is from February 24 to March 2, 
1979. Although the results for aerosol optical depth (Figure 3.2a) exhibit 
more fine structures than that of height field, their spatial and temporal 
variations are, as a whole, similar. 

The results for the 1980 data are given in Figure 3.3. In this case, the 

3 

coefficients a and b are found to be 22.3 and 1.81 x 10 , respectively. The 
value of the linear correlation coefficient is 0.525, derived from 853 data 
pairs. In contrast to the 1979 case, the linear correlation in 1980 is 

relatively poor. We have investigated this particular case in some detail 
and will discuss the results later. 

Figure 3.4 shows the results of the analysis for 1981. The number of 
data pairs used in this calculation is 647. The linear least-squares fitting 

3 

yields a constant 22.4 and a slope 2.47 x 10 . The calculated correlation 
coefficient is 0.753. Although this value is less than that of 1979 case, the 
linear correlation between aerosol optical depth and height is significant. 

The result of the analysis for January-March 1982 is displayed in Figure 

3.5. The number of data pairs employed in the calculation is 918. The values 

3 

for a and b are found to be 22.1 and 2.19 x 10 , respectively. The linear 
correlation coefficient is 0.812. Thus, there is a significant correlation 
between aerosol optical depth and 30 mb height. The degree of the linear 
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Figure 3.2. Projection plots of SAM II aerosol optical depth (a) and 
associated 30 mb height (b) as functions of longitude and 
time from January 24 - March 2, 1979. 
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Figure 3.3. Linear correlation between SAM II aerosol optical depth 
(above 20 km) and associated 30 mb height from January 
to March 1980. 
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Figure 3.4. Linear correlation between SAM II aerosol optical depth 

(above 20 km) and associated 30 mb height from January to 
March 1981. 
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Figure 3.5. Linear correlation between SAM II aerosol optical depth 
(above 20 km) and associated 30 mb height from January 
to March 1982. 
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correlation shown in this case is comparable with that of 1979. Table 3.1 
summarizes the results of the four described cases. 

To examine the January-March 1980 case in detail, the correlation 
analysis has been applied separately to the three individual months of 1980. 
The results are shown in Figures 3.6a to 3.6c and are summarized in Table 
3.2. As one may see, significant correlation values for the individual months 
are evident. A question is why the correlation is lower when it is evaluated 
with three months all together. It is found that the 30 mb heights for each 
of the three monthly data sets are substantially different from each other. 
This is apparent from Figure 3.7 where Figure 3.6 a-c are replotted on top of 
each other and are labelled differently for different months. The data pairs 
in January are shown by crosses (+), February by circles (o), and March by 
triangles( A) • The March data set as a whole is characterized by the highest 
30 mb height, the January data set by the lowest, while the February data set 
is intermediate. Since the atmosphere is very close to hydrostatic 
equilibrium, the greater 30 mb height implies generally the higher 
temperature. This feature is evident in Figure 3.8 which also shows the SAM 
II extinction and the corresponding temperature at 30 mb. The data pairs for 
different months are labeled as in Figure 3.7. Figure 3.8 indicates that the 
temperature in March as a whole is indeed the highest, January has the lowest 
temperature, and that February is intermediate. In conclusion, the 
correlation for January-March 1980 (Figure 3.3) is strongly dependent upon 
the temporal variations of the stratospheric temperature in the polar region, 
which do not seem to be as significant in 1979, 1981, and 1982. 

3.1.4 Remarks 

January to March aerosol extinction profiles from 1979 to 1982, derived 
from SAM II measurements and the associated meteorological information, have 


# 
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TABLE 3.1. Linear Correlation between SAM II Aerosol 
Optical Depth and Associated 30 mb Height 
from January to March. 


Year 

No. of data points 

a 1 

b* 

(x 10 3 ) 

0* 

1979 

1036 

22.1 

2.69 

0.834 

1980 

853 

22.3 

1.81 

0.525 

1981 

647 

22.4 

2.47 

0.753 

1982 

918 

22.1 

2.19 

0.812 


fa and b are the coefficients of a linear least-squares : 
fitting (y = a + bx) 

*0 is the linear correlation coefficient. 


* 
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1980 

JANUARY 

LINEAR CORRELATION, y = a + bx 
0 = 21.9, b= 2.05 x 10^ 

LINEAR CORRELATION COEFF = .748 



AEROSOL OPTICAL DEPTH. 1 .E-3 


Figure 3.6. Linear correlation between SAM II aerosol optical 
depth (above 20 km) and associated 30 mb height 
for 1980. (a) January. 
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1980 

FEBRUARY 

LINEAR CORRELATION, y = o + bx 

0=21.8, b= 2.97 x I0 3 
LINEAR CORRELATION COEFF = .760 



Figure 3.6 (continued). (b) February. 
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TABLE 3.2. Linear Correlation between SAM II Aerosol 
Optical Depth and Associated 30 mb Height 
(1980). 


Month 

No. of data points 

a* 

b* 

(x 10 3 ) 

0* 

January 

279 

21.9 

2.05 

0.748 

February 

279 

21.8 

2.97 

0.760 

March 

277 

22.8 

2.10 

0.813 

fa and b 
fitting 

are the coefficients of 
(y = a + bx) . 

a linear 

least-squares 
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Figure 3.7. Linear correlation between SAM II aerosol optical 
depth (above 20 km) and associated 30 mb height 
for 1980. Composite diagram. 
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been used to study the degree of linear correlation between aerosol optical 
depth (above 20 km) and the height of the 30 mb pressure level. Significant 
positive correlation is found for all the years except 1980. However, for 
1980, significant correlation does exist if the data is analyzed on a monthly 
basis. The different correlation behavior for 1980 when the three months are 
considered together can be attributed to the substantial temperature rise in 
the Arctic Region between January to March 1980. It should be noted that this 
is a period of stratospheric warming associated with amplifying planetary 
waves in the polar region. There is considerable variability in the 
manifestation of the warming events and amplifying planetary waves from year 
to year (Schorbel, 1978, Labit zke, 198 1 , Labitzke and Goretzki, 1982). The 
results in this study may indicate also that the change in the circulation 
pattern from Winter to Summer for 1979-1980 took place much later in 
comparison with what happened in the other three cases. This seems to be 
supported by the study of Labitzke (1981) and Labitzke and Goretski (1982). 
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4. TASK 3— TO INVESTIGATE THE INJECTION OF VOLCANIC MATERIALS INTO THE 
STRATOSPHERE, INCLUDING GLOBAL LOADING AND TRANSPORT (includes work 
carried out under Task Area 8) 

4.1 Volcanic Loading and Transport 

Several volcanic eruptions between 1979 and 1981 injected solid and 
gaseous material into the stratosphere, raising the mass of aerosol stored 
there. The resultant changes in 1 pm aerosol extinction were monitored by 
SAM II and SAGE I giving an almost global picture of the development and 
movement of the aerosol. The aerosol extinction was used to calculate 
approximate values for the aerosol mass loading in the stratosphere. From 
these values the global background mass loading has been determined for the 
prevolcanic period of 1979, together with an estimate of the seasonal 
variation of this quantity in each hemisphere. The changes in aerosol 
extinction in 1980 and 1981, following each volcanic eruption, have been used 
to calculate the mass of aerosol injected into the stratosphere by each 
eruption. The manner in which the material injected by these eruptions has 
divided itself between the Northern and Southern Hemispheres has permitted 
some general conclusions to be drawn about stratospheric circulation 
features . 

The analysis described above has been published in the Journal of 
Geophysical Research (Kent and McCormick, 1984), the abstract for this 
publication is given in Appendix G. Additional analysis of the same data set 
describing the rates of changes of aerosol mass loading following each 
eruption is given in Section 6.1 of this report. 
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5. TASK 4 — TO INVESTIGATE STRATOSPHERIC PLANETARY WAVES AND THEIR EFFECT ON 


TRANSPORT OF AEROSOLS, 0 3 and N0 2 

5.1 Aerosol Extinction and Temperature Changes Associated with Planetary 
Waves 

This section summarizes the analysis of behavior of the zonal mean 
aerosol extinction ratio and its relationship with the zonal mean temperature 
during the winter of 1978-1979 stratospheric warming. The detailed 
description of the investigation can be found in the report by Wang and 
McCormick (1985a, Appendix H). It is generally understood that the 
distribution of stratospheric aerosols is strongly influenced by aerosol 
microphysics and atmospheric dynamics. By examining the concurrent 
variations of aerosol extinction ratio and temperature, it is possible to 
recognize the key mechanism which controls the distribution of stratospheric 
aerosols. It is found that, based on SAM II aerosol extinction measurements, 
the zonal mean aerosol extinction ratio exhibits distinct variations during 
the January-February 1979 stratospheric sudden warming; and that the 
transport effect of planetary waves may have played a significant role in 
determining the distribution of values observed. 

5.2 Variations of Aerosol, 0 3 and N0 2 Concentrations during Stratospheric 
Warmings 

In the previous section, we have examined the effect of planetary waves 
on the distributions of stratospheric aerosols. In this section, the 
variations of stratospheric aerosols, and 0 3 and N0 2 concentrations during 
stratospheric warmings are investigated by utilizing both SAM II and SAGE I 
observations. Four specific data sets taken during the peak development of 
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planetary waves have been selected for this analysis* As mentioned earlier, 
both SAM II and SAGE I instruments provide about 15 measurements per day 

distributed along an approximately constant latitude at high latitudes. This 
measurement character istic allows us to examine the variations of the SAM II 
and SAGE I measured species that are associated with amplified planetary 
waves* It is found that there is a strong correlation between the 30 mb 
pressure height and the aerosol optical depth (at 1 year). Strong 

longitudinal gradients are found with the low values occurring within and 
wherever the vortex exists. Similar characteristics are also found to exist 

for the simultaneously observed SAGE 0^ and N0^ columnar density 

distributions. The detail of this analysis can be found in Wang and 
McCormick (1985b, Appendix I.) 

5.3 Ozone and Temperature Changes and Associated Fluxes During the 

Mid-February 1981 Stratospheric Warming 

5*3. 1 . Introduct ion 

In a sense, observation of the development of large-scale stratospheric 
ozone and temperature waves during the Winter at high latitudes would provide 
a unique opportunity for comparison with theoretical model prediction to 
assess current understanding of the behavior of the stratospheric ozone. The 
development of these large-scale waves as a response of the stratosphere to 
tropospheric disturbances and also the transport effects of these waves have 
been investigated using theoretical models (Hartmann and Garcia, 1979; Garcia 
and Hartman, 1980; Kawahira, 1982; Kurziga et al., 1984; Rood and Schoeberl, 
1983a, b, etc.). Meanwhile, studies based on satellite observations are 
becoming available (Barnett et al., 1975; Gille et al., 1980; Nagatani and 
Miller, 1984; Wu et al., 1985). In addition, ozone data based on the 
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Stratospheric Aerosol and Gas Experiment (SAGE) have been used in conjunction 
with the associated meteorological information to study the behavior of 
planetary waves of ozone and temperature, and their wave transports during 
the late February 1979 stratospheric warming (Wang, McCormick and Chu, 1983). 
Even though the analysis did not derive the detailed ozone budgets due to 
insufficient data, the results in Wang, McCormick and Chu, (1983) show many 
interesting features. More specifically, the ozone and temperature waves are 
shown to be in-phase in the lower stratosphere (below ^25 km), out-of-phase 
in the upper stratosphere (above ^38 km), with a transition region in 
between. In the upper stratosphere, an intense equatorward eddy ozone 
transport was shown to be accompanied by a poleward eddy heat transport 
during this specific warming. The analysis in Wang, McCormick and Chu (1983) 
also shows a rapid increase in ozone column density above 10 km over the data 
period analyzed. This increase was found to be primarily a response to the 
change of zonal mean ozone number density in the lower stratosphere. It 
should be understood that Wang, McCormick and Chu (1983) is a single case 
study. Since there is considerable variability in the manifestation of 

warming events and amplifying planetary waves from year to year (Schoeberl, 
1978; Labitzke, 1981; Labitzke and Goretzki, 1982), results from observations 
of different warmings are highly desirable. 

The purpose of this study is two-fold. First, the behavior of planetary 
waves of ozone and temperature and their wave transports during the middle 
February 1981 stratospheric warming will be studied in detail by utilizing 
the SAGE ozone measurements and associated meteorological data. Second, the 
analyzed results of this middle February 1981 warming will be compared with 
those of the late February 1979 warming described in Wang, McCormick and Chu 
(1983). 
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As will be shown, these two warming events correspond to two quite 
different stratospheric circulations at high latitudes. A comparison of the 
behavior of ozone and temperature waves between those two events is thus 
significant, and leads to enhancement of our understanding of their behavior 
during warming events. 

5.3.2 Data and approach 

Similar to Wang, McCormick and Chu (1983), the results of SAGE ozone 
measurements and the associated meteorological information provided by NOAA's 
National Meteorological Center are used in this analysis of the mid-February 
1981 warming event. In Wang, McCormick and Chu (1983), we have described the 
general characteristics of the SAGE ozone and the meteorological data sets 
relevant to the planetary wave studies. In this study, we will mention only 
specific aspects of these observations which are important to this particular 
analysis. As noted in Sec. 1.2, the AEM 2 satellite power system developed 
problems about 5 months after the SAGE instrument was launched. As a result, 
there were no sunrise measurements available after July 1979. The entire SAGE 
sunset observation period covers 34 months from February 1979 to November 
1981 (Figure 5.1). For the purpose of examining the behavior of winter 
planetary waves, it is desirable to use data obtained at high latitudes where 
the waves are most active. As shown in Figure 5.1, the SAGE measurement 
locations reach the highest latitudes on two occasions centered around 
January 1 and March 6, during the winter 1979-1980. Unfortunately, the 
highest latitude of the SAGE measurements [around January 1, 1980 is about 

46°N. As for the measurements near March 6, 1980], the highest latitude 

reached was 58°N, but data gaps are found in the retrieved profiles during 
this period. For these reasons, it was decided to examine the case of Winter 
1980-1981 instead. During this winter, there are 10 consecutive daily SAGE 
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ozone profiles from February 12 to 21, 1981 available for the analysis. The 
peak latitude near 55°N occurred on February 17, 1981. Table 5.1 indicates 

the' average latitude and the number of profiles of the SAGE daily 
measurements for these 10 days. The overall averaged latitude for this 
10-day period is about 52°N which is only slightly lower than the February 
1979 case discussed in Wang, McCormick and Chu ( 1983). It is important to 
note that this 10-day period covers approximately the second half of the 
major midwinter warming which occurred in February 1981 and reported by 
Labitzke (1982). 

In this analysis, the method for deriving the planetary waves of ozone 
and temperature and their fluxes using SAGE ozone and meteorological data 
sets is essentially the same as that described in Wang, McCormick and Chu 
(1983). For this reason, we will not repeat the discussion here. Since one 
of the purposes of this study is to compare the ozone and temperature waves 
and their eddy transports between the late February 1979 and the middle 
February 1981 warmings, a general description of the synoptical 
meteorological aspects of the northern high latitude stratosphere during 
these warming periods is included in the analysis. This will help to bring 
out some features of the atmospheric circulation and the thermal structure 
which are not obvious from an examination entirely in terms of their harmonic 
components. As will be seen in the next section, there are considerable 
differences in the stratospheric circulation pattern between these two 
warming events. 

5.3.3 Results and discussion 

In this section, we first show the development of the arctic circulation 
pattern for the middle February 1981 warming. For comparison, a description 
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TABLE 5.1. The Number of Profiles and the Averaged Latitude 
of SAGE Observations from February 12 to 21, 1981. 


Date 

(February) 

Number of 
Profiles 

Averaged Latitude 
(°N) 

12 

15 

49.47 

13 

14 

50.94 

14 

15 

52.05 

15 

15 

52.84 

16 

15 

53.34 

17 

15 

53.55 

18 

15 

53.50 

19 

15 

53.18 

20 

15 

52.59 

21 

14 

51.69 
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of the late February 1979 warming is included. The results of harmonic 
analyses of ozone mixing ratio, temperature and eddy meridional velocity, are 
given in the second part of this section, followed by a description of the 
horizontal meridional ozone and temperature transports by the planetary 
waves. In the last part of this section, the phase relationship between the 
waves of ozone mixing ratio, temperature and eddy meridional velocity is 
examined . 

5.3.3. 1 Synoptical meteorology of the high latitude stratosphere for the 

middle February 1981 warming 

Some aspects of the warming event of 1980-1981 have been discussed by 
Labitzke and Goretzki (1982). It began with the development of an Aleutian 
anticyclone which indicated the amplification of the wavenumber 1 disturbance 
in November 1980 — the so-called Canadian Warming. In December 1980, the wave 
activities were generally mild and the monthly mean 30 mb North Pole 
temperature was -82°C — the coldest temperature since the winter 1958-1959. 

This situation extended to the middle of January 1981. During this period, 
the high latitude Northern Hemisphere was characterized by the development of 
a vast quasi-circular polar vortex as shown in the 30 mb upper air map on 

January 14 1981 (Figure 5.2a). Thereafter, this polar vortex was disturbed by 
amplifying planetary waves, mainly wavenumber 1 (Figure 5.2b). By middle 

February, the high latitude zonal mean thermal structure had shown a reversal 
of the latitudinal temperature gradient, indicating the occurrence of a major 
midwinter warming. According to Labitzke and Goretzki (1982), the monthly 
mean 30 mb North Pole temperature had increased from -82°C in December 1980 
to -56°C in February 1981. Figure 5.2c shows the 30 mb upper air map on 

February 15, 1981, in which the remainder of the wavenumber 1 disturbance is 
still clearly evident on this particular day. It is interesting to note that 
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maps at 30 mb for (a) January 14, 198 





the wavenumber 1 height disturbance had shifted its phase by an amount close 
to 180° from January 25 to February 15. Note also that the cold center, which 
was located near the North Pole just before the development of planetary wave 
disturbances on January 14, 1981 (Figure 5.2a), had moved to 'V<60 O N latitude 
west of Canada by February 15, 1981. At the same time, warmer air masses had 
appeared in the polar region (Figure 5.2c). In the second half of February, 
the activity of planetary wave disturbances became much weaker and the polar 
vortex began to regain some of its strength (Figure 5. 2d). This is a stage 
corresponding to a period of the so-called "late winter cooling" 
(Labitzke and Goretzki, 1982). Later, there was a further development of 
wavenumber I disturbance which reached its peak intensity on about March 3, 
1981. Thereafter, the arctic stratosphere was gradually replaced by the 
summer circulation system. It is interesting to note that Labitzke ( 1982) 
had shown that this wavenumber I dominated winter disturbance took place 
during the westerly phase of the equatorial quasi-biennial oscillation (at 
50 mb) and not on the easterly phase as for most wavenumber 1 development. 
She further indicated that only wavenumber 1 winter events which occurred near 
sunspot maxima exhibit this unusual feature. For comparison purposes, it is 
essential to this investigation to also describe the circulation which 
occurred in the polar region during the late February 1979 warming. 

The synoptic aspects of the circulation of Winter 1978-1979 have been 
discussed by Noxon et al. (1979), and Syed and Harrison (1981) in their 
studies on the behavior of NC^ abundance during the warming of 
January-February 1979. The behavior of planetary wave disturbances of the 
temperature and height fields of this winter also have been investigated by 
Quiroz (1979) and Labitzke (1981) in detail. Before January 15, 1979, the 
polar stratospheric circulation system showed only mild fluctuations in these 
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fields except for a few days around December 8, J978. During these days, an 

Aleutian anticyclone developed indicating an amplification of height 

wavenumber 1. In the period January 15 to early March 1979, the northern 
polar stratospheric circulation system went through several interesting 
changes. First, there was the formation of an unusually strong Aleutian 
anticyclone which led to the amplification of height wavenumber 1 disturbance 
at 30 mb with the peak intensity occurring on about January 26. Around the 
same day, the stratosphere showed the first reversal of the meridional 
gradient of mean stratospheric temperature. Thereafter, this anticyclonic 
system weakened and the disturbed polar vortex tended to regain its strength 
with its center moving back toward the pole. Figure 5.3a shows the 30 mb 
arctic upper air map for February 9, 1979. This polar vortex became well 

established again by about 12 February. Beginning that day, the vortex 

exhibited an elongation of its circulation pattern approximately along the 
east-west axis, with the appearance of two low pressure systems (Figure 

5.3b). Meanwhile, two high pressure systems developed on each side of the 
elongated polar vortex (Figure 5.3b), indicating an amplification of height 
wavenumber 2. This development in the arctic stratosphere culminated in the 
splitting of the polar vortex (Figure 5.3c) and the reversals of the zonal 
mean arctic flow, as well as the meridional gradient of mean stratospheric 
temperature — a major final warming (Quiroz, 1979). By the end of February, 
these wavenumber 2 dominated disturbances began to decay. In the following 
month, there were only mild fluctuations of the height field in the form of 
wavenumber I. During the month of April 1979, the flow in the arctic winter 
stratosphere was gradually replaced by the summer circulation system. 

From the above discussion, it is clear that the development of the 
arctic winter disturbances in the two Winters 1978-1979 and 1980-1981 are 
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Figure 5.3. Constant pressure upper air maps at 30 mb for (a) February 9, 1979; (b) Feb 

1979; (c) February 25, 1979; and (d) March 2, 1979. isobaric altitude 

isotherms. Solid circles on (c) denote the SAGE measurement location 







rather different. In 1980-1981, the arctic disturbances involved primarily 
wavenumber 1 fluctuations during the entire winter season. In 1978-1979, 
however, the polar stratosphere showed disturbances involving both wavenumbers 
1 and 2 with each dominating during different periods. It should also be 
mentioned that in Wang, McCormick and Chu (1983), the period of the data set 
analyzed only covered the peak of the height wavenumber 2 development. On 
the other hand, the data set used in this analysis covers essentially the 

second half of the major midwinter warming which was during the decaying 
stage of a strong height wavenumber 1 disturbance. 

5. 3. 3. 2 Evolution of the planetary waves 

Figures 5.4a - 5.4c show the evolution of the temperature wavenumber 1 
and 2 disturbances and also the zonally averaged temperature during the 
mid-February 1981 stratospheric warming from February 12 to 21, 1981, 
respectively. In the case of the wavenumber 1 evolution (Figure 5.4a), the 
stratosphere exhibits three distinct layers. They are centered approximately 
at 22 km, 32 km, and 44 km, respectively. In the middle layer, a mild 
decrease in the wave amplitude is shown. The behavior of the waves in the 
lower and upper layers, is very similar, showing an intensification in the 
first 4 days, and a rather steady decay afterwards. The maxima of the wave 
amplitudes, approximately I2°C and 10°C in the upper and lower layers, 
respectively, appeared on February 15. The calculated amplitudes for 

temperature wavenumber 2 show generally a mild decay of the wave in the 
altitude region from 10 km to 50 km throughout the entire data period of this 
analysis (Figure 5.4b). From Figure 5.4a and 5.4b, it is clear that the 
temperature disturbances in this mid-February warming are dominated by 

wavenumber 1 activities. In the case of the zonal mean component (Figure 

5.4c), no significant variation is shown during this 10-day data period. 
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However, a 2°C increase in the mean temperature is noticeable between 28 km 
and 38 km in the period from the 5th day to the 10th day. 

In comparison with the case in Wang, McCormick and Chu (1983), the 
wavenumber 2 and mean temperature disturbances for this mid-February warming 
are smaller in magnitude. The wavenumber 1, however, is more intense 
overall, except in regions above 3 38 km where the two wavenumbers are 
comparable in magnitude. 

The derived amplitude of the first two harmonic components of the eddy 
meridional velocity are shown in Figures 5.5a and 5.5b, respectively. The 
wavenumber 1 shows an intensification from day 1 to day 3 above approximately 
20 km (Figure 5.5a). A peak value of 16.5 m/sec appeared on day 3 at about 37 
km. Thereafter, the wave decayed. Below about 20 km, the wave amplitudes 
are relatively small (Figure 5.5a). In Figure 5.5b, the result of wavenumber 
2 shows a decay of the wave throughout the data period of this study between 
10 km and 50 km. The strength of the wavenumber 1 disturbance (Figure 5.5a) 
is about 1.5 times that observed in Wang, McCormick and Chu (1983), while the 
wavenumber 2 strength of the disturbance is only about half. 

The first two harmonic components of ozone mixing ratio are given in 
Figures 5.6a and 5.6b. The noticeable features of the first component (Figure 
5.6a) are the appearance of a maximum at 44 km about February 15, 1981. The 
variation of ozone mixing ratio in a layer between 40 and 50 km seems to be 
associated with that of the temperature wavenumber 1 (Figure 5.4a). Below 
about 40 km, (Figure 5.6a), a relatively steady-state condition is observed. 
In the case of ozone wavenumber 2, the amplitude is only about half that of 
the wavenumber 1 (Figure 5.6b). Figure 5.6b also shows that most of the 
wavenumber 2 activities were in the upper stratosphere above about 30 km. 
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Figure 5.5. Evolution of the amplitudes (ms ) of meridional velocity 
waves during the middle February 1981 warming (^52°N); 
(a) wavenumber 1, (b) wavenumber 2; contour interval 
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Figure 5.6. Evolution of the amplitudes (ppmv) of the first two ozone 
waves, and the zonal mean ozone during the middle February 
1981 warming («^52°N); (a) wavenumber 1; (b) wavenumber 2, 
contour interval 0.2 ppmv, and (c) zonal mean ozone, contour 
0.8 ppmv. 




Below 30 km, the wavenumber 2 is relatively weak. The evolution of zonal 
mean ozone mixing ratio is given in Figure 5.6. Below about 30 km (Figure 
5.6c) little change is observed in the mean ozone mixing ratio during the 
data period of this study. In the layer between 30 km and 42 km, there is a 
slow increase in the ratio in the late half period. Above 42 km, there is a 
noticeable decrease in the ozone mixing ratio in the entire data period of 
this study. In comparison, the data shown in Figures. 5.6a and 5.6c are, in 
many respects, similar to those shown in Wang, McCormick and Chu (1983, 
Figure 5.6.) In both cases, the ozone wavenumber 1 disturbances are observed 
to be more intense in the upper stratosphere component. Slight increases are 
also observed in the zonal mean ozone mixing ratio lower stratosphere below 
30 km. Particular noteworthy is the fact that ozone wavenumber 1 
disturbances in the upper stratosphere are clearly correlated to that of 
temperature wavenumber 1 in these two different warmings. This correlation 
between temperature and ozone wavenumber 2 amplitudes, however, is less clear 
in the analysis of these two warmings. 

In order to examine the changes in the columnar ozone above 10 km during 
this mid-February 1981 stratospheric warming, the evolution of zonal mean 
ozone number density in this event is also derived based on the SAGE ozone 
data set. Figure 5.7a shows the computed results. From Figure 5.7a, it is 
evident that there is an increase in ozone concentration in the lower 
stratosphere centered at altitude 20 km. Above 24 km (Figure 5.7a) no 
apparent changes in the mean number density are observed. It is interesting 
to note that the behavior of zonal mean ozone density in this mid-February 
1981 warming is very similar to the late February 1979 event (Wang, McCormick 
and Chu, 1983). The associated time variation in the ozone columnar density 
is given in Figure 5.7b. It shows a monotonical increase in value except for 
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Figure 5.7. Evolution of (a) zonal mean ozone number density, contour 
interval 0.8, scale by 10~ 12 in unit cnT 3 ; and (b) ozone 
columnar density at altitude 10 km during the middle 
February 1981 warming (^52°N). 
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the first two days. This increase in ozone columnar abundance is primarily a 
response to the increase in the ozone concentration in the lower stratosphere 
as illustrated in Figure 5.7a. A similar situation was also found in the late 
February 1979 stratospheric warming event (Wang, McCormick and Chu, 1983). In 
comparing the ozone columnar density of this analysis with that described in 
Wang, McCormick and Chu (1983), it is found that the overall values in the 
late February 1979 event are generally greater than those occurring in this 
mid-February 1981 warming. 

5. 3. 3. 3 Horizontal ozone and temperature transports by planetary waves 

The computed eddy ozone transports associated with the first two wave 
components in the mid-February 1981 warming are displayed in Figures 5.8a and 
5.8b. The remarkable features of the wavenumber 1 transport (Figure 5.8a) are 
the clearly defined altitude regions for the poleward and the equatorward 
ozone transports throughout the entire data period of this analysis. Above 
40 km, a distinct layer of equatorward eddy ozone transport is observed 
(Figure 5.8a), with a maximum value of 4.5 ppmv ms * occurring on February 
14, 1981 at 44 km. Below about 28 km, an equatorward transport is also 
observed (Figure 5.8a), but it is much weaker in comparison with the one 
above 40 km. In the layer between 28 and 40 km, the eddy ozone transport is 
mainly toward the North Pole. A slight increase in the transport is observed 
in this region during the first 4 days which decayed thereafter. A maximum 
poleward transport of VJ ppmv ms * appeared at 36 km on the same day 
(February 15, 1981) that the equatorward transport reached its peak value in 
the upper stratosphere. 

An analysis of the ozone transport associated with wavenumber 2, 
indicates a poleward transport below 40 km and equatorward transport above 40 
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Figure 5.8, Time variations of eddy ozone flux (ppm ms“^). Due to 

(a) wavenumber 1, contour interval 1 ppm ms“l , and 

(b) wavenumber 2, contour interval 0,5 ppm ms“l. The sum 
of the first three waves is given in (c) , contour interval 
1.5 ppm ms~l. (+ poleward, - equatorward, ) 
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km in the first 2 to 3 days (Figure 5.8b). Later, the equatorward transport 
system in the upper stratosphere is observed to increase slightly reaching a 
peak value of 3 ppmv ms * at 46 km on February 15, 1981. On February 14, 1984 
(Figure 5.8b), a change in transport direction is observed within the layer 
between approximately 24 km and 38 km from poleward to equatorward. After 
February 15, the wavenumber 2 eddy ozone transport in the entire layer is 
observed to decay. From Figures 5.8a and 5.8b, it is evident that the 
wavenumber 2 eddy ozone transports are generally weaker than wavenumber 1 in 
this particular warming. The sum of the eddy ozone transport due to the 
first three wave components is given in Figure 5.8c. Since the transport due 
to wavenumber 1 is relatively strong, the transport shown in Figure 5.8c 
exhibits features similar to those shown in Figure 5.8a. By comparing Figure 
5.8c with the corresponding results presented in Wang, McCormick and Chu 
(1983), one will notice that they exhibit generally similar features. In 
both cases, the ozone mixing ratio transports exhibit three distinct layers, 
namely, two equatorward transport layers, one in the upper stratosphere, and 
one in the lower stratosphere, and a poleward transport layer in the middle 
stratosphere. There are, however, slight differences in the vertical 
extension and the center altitude of these layers for the two warming cases. 

Since eddy ozone mass transport is directly related to changes in total 

ozone amount, the contribution of the first two wave components to this 

transport during the mid-February 1981 warming are analyzed (Figures 5.9a and 

5.9b). Below about 26 km, the development of an intense equatorward 

wavenumber 1 transport is observed (Figure 5.9a). This transport reached a 

12 -3 -1 

peak intensity of 2.7 x 10 molecules cm ms at 16 km about February 17, 
1981. Above 26 km, the wavenumber I transport is relatively weaker than that 
below. Nevertheless, a layer of poleward mass transport between 26 km and 38 
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Figure 5.9. Time variations of eddy ozone mass transport (molecule cm ms ). 

Due to (a) wavenumber 1; and (b) wavenumber 2, contour interval 
0.25 molecule cnT^ms”!. The sum of the first three waves 
is given in (c) contour interval 1 molecule cm~3. Scaled 
by 10~12 (+ poleward, - equatorward. ) 
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km and also a layer of equatorward transport above 38 km are noticeable 

(Figure 5.9a). The calculated results of wavenumber 2 ozone mass transport 

indicates a contribution that is generally much smaller than that of 

wavenumber 1 except in regions below 22 km during the first 2 to 3 days 

(Figure 5.9b). In these regions, poleward eddy transport is observed during 

those days. The sum of the first 3 components of eddy ozone mass transports 

is displayed in Figure 5.9c. It is obvious that an intense equatorward 

transport developed below approximately 26 km. This development seems to be 

associated with the decaying of a poleward eddy ozone mass transport in the 

layer between 26 km and 40 km as revealed from the first 3 days of this data 

period. Although this poleward transport became rather weak by the third 

day, a thin layer of poleward transport centered at altitude 34 km is still 

noticeable in the rest of the data period. From Figure 5.9, it is clear that 

the wavenumber 1 component plays a dominant role in the eddy ozone mass 

transport in this mid-February 1981 warming. As indicated in Figure 5.9c, 

the peak of the equatorward transport occurred around February 16, 1981 at 

~16 km. The vertical integrated mass flux of ozone (10-30 km altitude) on 

20 -2 

this day is found to be (3.2 0.08) x 10 (molecules/cm )(cm/sec) to the 

south. 

In Wang, McCormick and Chu (1983), the ozone and temperature eddy 
transports in the upper stratosphere (above 35 km) were shown acting in 
opposite directions in both wavenumber 1 and 2 components. Of particular 
interest is that both the wavenumber 1 ozone and temperature eddy transports 
show a change in their direction on about the same day. To see whether such 
a feature exists in the middle February 1981 warming, the eddy transports of 
temperature have been calculated. The results are shown in Figure 5.10. The 
most distinct feature of the temperature wavenumber 1 transport (Figure 
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5.10a) is the development of a poleward transport above about 30 km with the 
peak at about 44 km. This peak reached its maximum intensity (^67 km/sec) on 
approximately February 14, and decreased steadily thereafter. Below about 30 
km, only mild equatorward eddy heat transports are observed (Figure 5.10a). 
In the case of wavenumber 2, a decrease of the poleward heat transport is 
generally observed throughout most of the altitude range of this study during 
this middle February 1981 warming (Figure 5.10b). The net eddy heat 
transports (shown in Figure 5.10c) exhibit a rather similar time variation as 
that of wavenumber 1, with a distinct layer of poleward eddy heat transport 
above 30 km centered at about 40 km, particularly on February 14. It should 
be noted that the daily values of the temperature difference between 80°N and 
50°N at 30 mb pressure level exhibit the second peak of the temperature 
reversal on about February 14, 1980 (Labitzke, 1982). This peak occurs about 
the same time at which the eddy heat transports began to change its direction 
from poleward to equatorward. As mentioned earlier, the wavenumber 1 ozone 
transport during the period of the data set in this study shows the peak of 
equatorward transport at an altitude of 44 km. In addition, this peak 
reached its maximum intensity on February 14. Thus, the observed wavenumber 1 
of ozone and temperature transports are negatively correlated during this 
mid-February 1981 warming. In the case of wavenumber 2, the behavior of 
ozone and temperature transports showed similar negative correlation above 40 
km. In particular, in this altitude region, they both changed their 
transport directions sometime between February 20 and 21. 

5. 3. 3.4 Phase relationships between the eddy fields 

The phase relationship between ozone, temperature, and meridional 
velocity waves of wavenumber 1 component on February 13, 1981, is shown in 
Figure 5.11a. The horizontal bars, as in Wang, McCormick and Chu (1983), are 
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Figure 5.11. Phase relationship between ozone mixing ratio (solid line), temperature (dashed line), and 
eddy meridional velocity (solid and dashed line) waves during the middle February 1981 
warming ('-~52°N). Phase increases westward: (a) wavenumber 1, February 13, 1981; 

(b) wavenumber 2, February 13, 1981; (c) wavenumber 1, February 17, 1981; and 
(d) wavenumber 2, February 17, 1981. 



the computed uncertainty of the calculated phase based on the given 
uncertainty of the SAGE and meteorological data. On the one hand, Figure 
5.11a shows that the ozone and temperature waves are nearly out-of-phase 
above approximately 42 km. Below about 27 km, on the other hand, they are 
closely in-phase. In the region between approximately 27 km and 42 km, a 
change from the in-phase relationship below 27 km to an out-of-phase one 
about 42 km takes place. As to the phase relationship between ozone and 
meridional waves, it is observed that they are generally out-of-phase below 
about 25 km and above approximately 43 km (Figure 5.11a). They become closely 
in-phase at altitudes centered at 30 km. This phase relationship explains 
the vertical variation of the wavenumber 1 ozone transport which was observed 
on February 13, 1981 (Figure 5.8a). The phase relationship between 

temperature and meridional waves can also be seen in Figure 5.10a. In the 

region above 32 km, their phase difference is less than tt/2, and is greater 
than tt/2 in the region below. This feature explains the two distinct 

regions of different heat transport directions which was observed on the 

corresponding day (Figure 5.10a). The phase relationship of wavenumber 2 
components on February 13, 1981 is shown in Figure 5.11b. The ozone and 

temperature waves are observed to be out-of-phase above 42 km. Below 38 km, 
they are generally in-phase. A similar phase relationship also can be found 
to exist between ozone and meridional velocity waves. The wavenumber 2 
component of the temperature and meridional velocity waves are observed to be 
closely in-phase throughout the entire altitude range of this study on 

February 13, 1981. This phase relationship between temperature and meridional 
waves accounts for the poleward heat transport at all altitudes on this 

particular date. The phase relationship between ozone, temperature, and 

meridional velocity waves on February 17, 1981 are shown in Figures 5.11c and 
5. lid, for wavenumber 1 and wavenumber 2, respectively. As can be seen, the 
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phase relationships are generally similar to those on February 13, between 
the corresponding wave components. 

The detailed time variation of the wavenumber 1 phase relationship over 
the period of the data set at altitudes 44 km and 26 km are shown in Figures 
5.12a and 5.12b, respectively. The former shows the typical phase 
relationship in the upper stratosphere and the latter depicts that in the 
lower stratosphere. It is evident that an out-of-phase relationship between 
ozone and temperature waves in the upper stratosphere and a nearly in-phase 
relationship in the lower stratosphere occurred throughout during this 
mid-February 1981 warming. Similar phase relationships can also be found for 
the wavenumber 2 components (Figures 5.12c and 5.1 2d). The evolutions of the 
phase of meridional velocity waves are also shown in Figures 5.12a to 5.12d 
for comparison. 

5.3.4 Summary 

We have analyzed the planetary waves of ozone mixing ratio, temperature and 
meridional velocity near 52°N during the middle February 1981 stratospheric 
warming by utilizing SAGE ozone data and the associated meteorological 
information. In addition, the eddy ozone and temperature fluxes during this 
warming were also examined. The result large temperature wavenumber 1 
developments in the upper stratosphere centered at altitude 44 km and also in 
the lower stratosphere with a peak at about 22 km. The maximum intensities 
of about 12°C and 10°C in the upper and lower stratosphere, respectively, 
occurred on approximately February 15, 1981. Comparatively, the corresponding 
wavenumber 2 component is much weaker in intensity, especially in the upper 
stratosphere. Nevertheless, the wavenumber 2 activities in the lower 
stratosphere, centered at 22 km, are clearly visible in the computed results 
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Figure 5.12. Time variation of the phase relationship between ozone mixing 

ratio (solid line), temperature (dashed line), and eddy meridional 
velocity (solid and dashed line) waves during the middle 
February 1981 warming (^ 52°N) . Phase increases westward: 

(a) wavenumber 1 , 44 km altitude; (b) wavenumber 1, 26 km; 

(c) wavenumber 2, 44 km; (d) wavenumber 2, 26 km. 




and were observed to decrease during almost the entire data period of this 
study. With respect to the meridional velocity waves, the peak intensity of 
the wavenumber 1 component occurred at an altitude about 37 km on February 

14, 1981. In contrast, the corresponding wavenumber 2 component exhibited an 

overall decline of the wave amplitude during this middle February 1981 
warming. The analysis of the ozone mixing ratio wave during this period 
exhibits a distinct development of wavenumber 1 in the upper stratosphere 
centered at 44 km with the peak intensity occurring on February 15, 198 1. In 

the lower stratosphere, the ozone wavenumber 1 was much less active and 
remained rather steady throughout the entire period of this middle February 
warming. This trend also seems to be the case for the wavenumber 2 ozone 

mixing ratio waves. With respect to ozone eddy transports, the results show 
a layer with poleward transport between approximately 30 km and 40 km. This 
transport was found to be decreasing during this warming period. Above 40 
km, ozone was being transported equatorward with the peak intensity occurring 
on February 15 at 44 km. The results also indicate that the wavenumber 1 
activity is the dominant component in this ozone transport. In the case of 
eddy temperature transport, the results show large poleward fluxes above 
about 30 km. The peak of the fluxes occurred on February 14, 1981 at 

altitude 44 km. Similar to the case of ozone transport, wavenumber 1 
component is found to be primarily responsible for the temperature 

transport . 

Furthermore, the results of this warming have also been compared with 
those of the late February 1979 warming. For completeness, the synoptic 
aspects of the large scale disturbances during the Winter 1980-1981 were 
discussed. Similar discussions on that of the Winter 1978-1979 were also 
included for comparison purpose. It is interesting to note that, in contrast 
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to the late February 1979 warming which involves disturbances with an 
outstanding wavenumber 2 component, this middle February 1979 warming is 
associated primarily with large amplifying wavenumber 1 disturbances (Figures 
5.2c and 5.3c). Despite this distinct difference in the planetary wave 
activities between these two warming events, many common features are clearly 
evident. First of all, the results show that the development of ozone 
wavenumber 1 component in the upper stratosphere was closely associated with 
that of corresponding temperature waves. More specifically, they are 
negatively correlated. This feature is also noticeable in the case of 

wavenumber 2 component although it is not as predominant as wavenumber 1. 
Furthermore, in both warmings, there was rapid increase in the 

zonally-averaged ozone column abundance. This increase was shown to be a 
response to the increases of the zonal mean ozone number density in the lower 
stratosphere centered at altitude 20 km. With respect to the ozone mixing 
ratio transports associated with the waves, the results generally show the 
equatorward transports in the upper stratosphere and poleward transports in 
the middle stratosphere. In the lower stratosphere, the transport is much 
less intense. The equatorward ozone transport in the upper stratosphere is 
found to be accompanied by poleward heat transport during both warming 
events. As to the phase relationship between ozone and temperature waves, 
the computed results of this warming event also show three distinct layers, 
i.e,, an out-of-phase relationship in the upper stratosphere and an in-phase 
relationship in the lower stratosphere with a transition layer in between. 

Finally, in contrast to the warming case studied by Wang, McCormick and 
Chu (1983), which involves an outstanding wavenumber 2 activity, this 
analysis was devoted to a warming dominated by wavenumber 1 development. 
Although the data period covers only 10 days, it includes the second half 


98 



part of the major midwinter warming according to Labitzke and Goretzki 
(1982). Furthermore, like the analysis of Wang, McCormick and Chu (1983), the 
results of this investigation can be useful for modeling the winter 
circulation in the northern high latitude stratosphere and for simulating the 
planetary wave development and their transport effects. 
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6. TASK 5 — TO INVESTIGATE THE FORMATION AND EVOLUTION OF AEROSOLS INCLUDING 


THE STUDY OF PHYSICAL AND CHEMICAL PROCESSES 

6. 1 Post-Volcanic Decay Rates 

6.1.1 Introduction 

Stratospheric aerosol concentrations are strongly dependent upon 
volcanic injection of sulfur dioxide followed by chemical conversion into 
sulphate and growth of sulfuric acid aerosols. Theoretical simulations and 
experimental measurements of aerosol optical properties show two distinct 
phases (Deirmendjian 1973; Deepak 1982). There is an initial growth phase 
lasting 3-6 months, followed by a slower quasi-exponential decay over a 
somewhat longer time period. Comparison of theoretical models and direct 
measurements have suffered in the past from the fact that the experimental 
observations are usually made from a single location on the ground or, if 
from an airborne platform, are limited to a relatively short time period. In 
both cases, the experimental results are strongly affected by poorly known 
aerosol movements and may not represent well the global variation in 
stratospheric aerosol concentration. 

The advent of the SAM II and SAGE satellite systems with global 
capability for measuring stratospheric aerosol extinction (McCormick et al. 
1979) has enabled a much more accurate picture to be obtained of the aerosol 
mass variation following volcanic injection (Kent 1982; Kent and McCormick 
1984). It is thus of interest to make a more detailed comparison of these new 
experimental results with the predictions of the theoretical models and to 
use this comparison to determine the principal contents of the aerosol growth 
and decay. The approach used in this paper has been to search for a simple 
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empirical function that provides a good fit to an available theoretical model 
(Turco et al. 1982). This function, containing adjustable constants has then 
been fitted to the global optical depth time sequence as obtained from the 
SAGE/SAM II measurements between February 1979 and November 1981. The use of 
a fitted function, rather than an apparently simpler direct measurement from 
the optical depth sequence, is made necessary by the fact that five eruptions 
of significance occurred during this time period with overlapping and 
superimposed stratospheric effects. 

6.1.2 Derivation of empirical function and fitting to theoretical curves 

A recent theoretical study, showing the predicted long-term variation in 
stratospheric aerosol optical depth following volcanic injection, is that of 
Turco et al., 1982. This study simulates the physical and chemical processes 
associated with the stratospheric aerosol following the eruption of Mt. St. 
Helens in May 1980. The resultant optical depth variations (at wavelength of 

0.55 pm) were determined as a function of time after the eruption for 
different injected quantities of S0 2 and H 2 0. These variations are shown by 
the solid lines in Figure 6.1. An empirical analytic function, providing a 
good fit to these theoretical curves, was searched for subject to the 
following criteria: 

1. It should not contain more than two or three floating parameters. 

2. It should fit both the rising and falling portions of the growth 
curves . 
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3, It should be subject to easy computation and give unambiguous figures 
for the time between the eruption and the maximum optical depth and the initial 
rate of decay following the maximum. 

The functional form adopted was 

f(t) - A + B t n e -t/t ° 


( 1 ) 


where 


t = time after the eruption 
f(t) ■* optical depth 

A * a constant (background level prior to eruption) 

and B, n and t Q are empirical constants to be determined from the fitting. 

Examination of this equation shows that the time between the eruption 
and maximum optical depth is given by 


t 

m 


n t 

o 


( 2 ) 


and that the maximum optical depth 

f - A + B t n e _t n/ t o 
m m 

( 3 ) 

The time, t., taken for the increase in optical depth to decay by a factor 
d 

]/e from its maximum value is given by the equation 
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1 + I/n 


(4) 


(t./t ) - £n(t./t ) - 

a ID O m 

For large t, the decay becomes quasi-exponential with time constant t Q . 

Equation ( 1) has been fitted to both the theoretical and experimental 
data using the Levenberg-Marquardt algorithm (Brown and Dennis, 1972) for a 
non-linear least squares approximation, the calculations being carried out on 
the Control Data computer system at NASA/LaRC. In the case of fitting to the 
theoretical model, Turco et al. f s theoretical curves have been read at 
0.75-month intervals from 0 to 9 months, 1.5-month intervals from 9 to 12 
months and 3-month intervals from 12 to 27 months. The resultant curves for 
the fitted empirical function are shown by the dashed lines in Figure 6.1, 
the derived parameters being listed in Table 6.1. From inspection of Figure 
6.1, it can be seen that the theoretical models and the fitted analytic 
functions agree very well. Table 6.1(b) shows the difference between the 
three main parameters, f , t and t,, derived from the fitted functions and 
directly measured from the theoretical curves. With the exception of the 
time delays for Model (1), where there is a three-week error in the 
determination of time of maximum, the errors are extremely low. Taking into 
consideration that the purpose of using an analytical model function is to 
derive these parameters from experimental data which itself will contain 
considerable uncertainties, these differences are considered acceptable. It 
should be noticed that the empirical formula has been chosen to fit the first 
two years variation after an eruption and is not necessarily representative 
of the behavior over a longer time period. 

6. 1.3 Fitting to experimental data 

The experimental data set used is that obtained by the SAGE and SAM II 

satellite systems between February 1979 and November 1981. The data set has 

been used by Rent and McCormick (1984), to derive the global mean 
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TIME (MONTHS) 

Figure 6.1. Modeled variations in stratospheric aerosol layer optical 

depth (at 0.55 pm) following volcanic eruptions with different 
amounts of injected material (Turco et al., 1982). Solid 
curves show results of theoretical calculations, dashed curves 
the fitted empirical functions. 
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TABLE 6.1(a). Constants Derived from Empirical Fitting to Theoretical Models in Figure 
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stratospheric aerosol optical depth at 1.0 jim measured from 2 km above the 
tropopause for the Northern and Southern Hemispheres over this time period 
and is fully described by them. Large variations in optical depth were 
observed which were attributed to one or the other of five volcanic eruptions 
which affected the stratosphere during this period. Figure 6.2 shows the 
mean global optical depth as a function of time and also shows the times of 
eruption of these five significant eruptions. Prior to the first eruption, 
that of Sierra Negra in November 1979, the mean global optical depth is 
almost constant. Following each eruption, the optical depth increases to a 
maximum after 4-5 months, following which a decrease is observed. The 
eruptions have occurred sufficiently close together in time that the 
stratospheric effects overlap. This is particularly true for Alaid and Pagan 

which we shall not be able to separate analytically and for St. Helens and 
Ulawun which are separable. In two places, just after the eruptions of St. 
Helens and Ulawun, double values are shown for the optical depth. One of 
these corresponds to the total global average. The other has been obtained 
by averaging over a hemisphere (extending from the North to South Poles), 
which, because of the still localized nature of the stratospheric injection 
from the recent eruption, is unaffected by that eruption. Error bars 
representing the uncertainty in the mean optical depth are shown in the 
figure; the sources of error have been discussed in detail by Kent and 
McCormick (1984). 

The maximum optical depth, the time between each volcanic eruption and 
that of the corresponding maximum optical depth and the decay time for the 
injected aerosol has been derived for each eruption by fitting a function of 
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TABLE 6.2. Key Background and Volcanic Parameters Determined from Empirical Fitting. 
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the form 


i = 5 

f(t) = A + I B. (t - t.)"i e 
L 1 1 

i = 1 


(t tj.) / ^oi 


(5) 


to the data points shown in Figure 6.1. In this equation, A is constant 
representing the prevolcanic background condition and the suffix i refers in 
sequence to each of the five volcanic eruptions. t^ refers to the time of 
the ith eruption and n^ and t^. have the meanings defined earlier in relation 
to the ith eruption. It is, in addition, assumed that the contribution to 
the total optical depth of any volcano, prior to its eruption, is zero. 
Because of the impossibility of separating the effects of Alaid and Pagan 

within a global average, it has been necessary to set = B^, n^ = n^ and 

t . = t _ in order to stabilize the solution. Equation 5 thus contains 13 

unknown parameters representing the background and the various eruptions. 
Despite this, a unique solution was readily obtained which has been used to 
derive the peak optical depth and appropriate time constants for each 
eruption. Data for the peak optical depths and their use to estimate the 

mass loading for each volcano has been presented by Kent and McCormick (1984) 

and will not be discussed in any detail here. 

Although use of Equation (5) permits the derivation of a unique set of 
parameters, it is not obvious how the errors in the values for optical depths 

shown in Figure 6.2 should be used to derive corresponding errors in the 

derived parameters. This problem was handled by a statistical approach. The 
optical depth errors were used, in conjunction with a random number 

generator, to randomly increase or decrease each measured optical depth value 
from its mean value by an amount equal to its standard error. Each set of 

new optical depth values thus obtained was used in conjunction with Equation 
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(5) to determine a set of parameters for each volcanic eruption. This 
process was repeated a large number of times (approximately two hundred) 

using different sets of random numbers and the mean and standard deviation 
found for the total set of values for each parameter. These are shown in 
Table 6.2. The best estimate of f(t) [Equation (5) is shown by the solid line 
in Figure 6.2.] 

The fit between the analytic curve and the experimental data shown in 
Figure 6.2 is seen to be extremely good, the curve lying within the limits of 
the error bars for most data points. This fit indicates that there is a 
close relationship between the theoretical and experimental results and that 
both are well simulated by Equation (1). In noting this fact, it should be 
remembered that the theoretical optical depth curves were calculated for a 
wavelength of 0.55 pm whereas the experimental measurements were made at a 
wavelength of 1.0 pm. The two are not equivalent and their relationship 

depends on the aerosol size distribution and composition of the volcanic 

aerosol, which will certainly change during the months following an 
eruption. Thus, even though the analytic forms used for the growth and decay 
curves are the same, the parameters of the curves might be expected to differ 
somewhat in the two cases. 

Values for t and t^ in Table 6.2 may be compared with each other and 
also with the theoretical values, in Table 6.1. The errors in Table 6.2 may be 
compared with those shown in Table 6.1. The main conclusions to be shown are: 

1. With the exception of the Ulawun eruption, whose effects are 

difficult to separate from those of St. Helens, errors in the calculated 
parameters are less than about 1 0%. 

2. The time between the eruption of each volcano and its maximum effect 

on the stratosphere, t is between 3-1/2 and 4-1/2 months for all the 

m 

eruptions studied, agreeing quite well with the values in Table 6.1(a). 


Ill 


3. The 1/e decay time, t. shows more variability and more error than t 

a m 

with its average value being about 6 months. This is somewhat smaller than 
the values in Table 6.1(a). 

4. Errors in f , t and t, in Table 6.2 are, in general, somewhat larger 

m m d 

than those in .Table 6.1(b) and are, therefore, not greatly effected by the 
inaccuracies in the empirical model. 

6.1.4 Comparison with other observations and conclusions 

Many observations have been made of the remote optical and in situ 
changes produced in the stratosphere by volcanic eruptions. These 

measurements and the subsequent estimates of the aerosol mass loading and its 
rate of change suffer from the disadvantage of having only limited sampling 
in either space or time or both. This can effect very strongly the observed 
time delay, t to maximum stratospheric change, as the maximum effect 
observed from a single point on the ground may be produced by the passing of 
a single intense cloud of material which is unrepresentative of the global 
distribution. To a lesser extent, the decay rate is also effected, being 
strongly dependent upon the closeness of the observing point to any 
stratospheric aerosol reservoirs and to the seasonal movements within and 
between the hemispheres. Data comparisons may thus be regarded as 
semi-quantitative rather than exact. While general conclusions may be drawn 
regarding the behavior of the stratosphere following an eruption, any exact 
study, such as distinguishing the characteristics of one eruption from 
another, must be regarded as difficult. 
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Table 6.3 lists some of the observations made of stratospheric effects 
following volcanic injection, together with estimates of the relevant time 
constants. The list is not exhaustive, but includes both observations made 
on earlier eruptions and observations made by other workers on the eruptions 
discussed in the present paper. Latitudes at which the observations were 

made, as well as those of the eruptions, are shown in the table. In general, 
quantitative terms, the values for the two major time constants, are 
consistent with each other and with the values shown in Table 6.2. In detail, 

the agreement is not good and illustrates the comments made above concerning 

the problems of determining volcanic time constants from a single or a small 
number of ground observation points. Thus, the low values of 1-1/2 to 2 

months, observed for the time delay to maximum effect at some stations, are 

almost certainly due to the passage of single intense volcanic clouds over 

those stations. The range of values observed by Hofmann and Rosen (1983, 

1984) for El Chichon arises from the different height ranges and different 

parameters measured. This list illustrates the difficulty of extrapolating 
from local to global terms. 

The main conclusion to be drawn from the analysis presented in this 
paper is that it is now possible, given the global coverage of satellites, 
such as SAM II and SAGE, to establish fairly precisely the rates of change in 
the stratospheric aerosol concentration following a significant volcanic 
eruption. The empirical model used here represents a useful technique for 
doing this particularly when several eruptions occur within a limited time 

frame, such that their atmospheric effects overlap. 
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7. TASK 6— TO INVESTIGATE THE DETERMINATION OF STRATOSPHERIC MOLECULAR AND 
AEROSOL OPTICAL PROPERTIES USING MULTIWAVELENGTH EXTINCTION AND 
BACKSCATTER DATA 

7.1 Retrieval of Aerosol Properties from Aerosol Extinction Ratios at 0.45 pm 
and 1 . 0 pm 

The possibility of retrieving the size distribution of aerosol particles 
in the stratosphere by measuring the extinction of solar radiation traversing 
the aerosol medium at wavelengths 0.45 urn and 1.0 pm has been explored in 
this study. The aerosol size distributions used consisted of nine analytical 
expressions, including most stratospheric aerosol models used by 
investigators. The aerosol compositions used were supercooled sulfuric acid 
droplets, with different weight percentage of I^SO^ in the aerosol. 

In each aerosol size expression used, there are two unknowns, namely, 

the constant N, which determines the number concentration of aerosol 

particles and the variable parameter, r , which determines the mode radius of 

S 

the aerosol size distribution. To eliminate one unknown, we calculated R, 
the ratio of the volume extinction coefficient of solar radiation by aerosol 
particles at wavelength X - 0.45 pm to the volume extinction coefficient at 
wavelength X ■ 1.0 pm. Results of this parametric study showed that R was 
not very sensitive to either the composition or the radii limits of 
stratospheric aerosol under consideration, but was quite sensitive to the 
value of the variable parameter governing the mean radius. In addition, 
there was a one-to-one correspondence between the extinction ratio R and the 
variable parameter in all nine aerosol size distributions. Consequently, the 
experimental value for R could be utilized to retrieve the variable parameter 
from the plotted curves generated in this study. After the variable 
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parameter had been retrieved, the expected aerosol extinction coefficient for 
N*1 could be calculated using Mie scattering theory and the true value of N 
for aerosol particles in the atmosphere was then given by the ratio of the 
measured aerosol extinction coefficient to the expected aerosol extinction 
coefficient for N-l. 

This work has been published in Applied Optics (Yue and Deepak, 1983). 
The abstract is included in Appendix J. 

7.2 Multiwavelength Aerosol Extinction Models 

Between 1979 and 1982, there were several volcanic eruptions that 
injected material into the stratosphere. This series culminated in the very 
large eruption of El Chichon whose effects are only now disappearing (in 
1986). In view of the great changes to the stratospheric aerosol layer that 
were observed during this period, it is of interest to ask in what manner 
these changes have affected the multiwavelength extinction characteristics of 
the aerosol. A preliminary study was made of these changes, using the seven 
wavelengths listed in Table 7. I (SAGE II wavelengths) and three stratospheric 
models. The models are as follows: 

MOD 4: A prevolcanic model (based on measurements by Gras and Michael, 
1979) 

MOD 5: A post Mt. St. Helens model (based on measurements by Hofmann 
and Rosen, 1982) 

MOD 6: A post El Chichon model (based on measurements by Hofmann and 
Rosen, 1984) 
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TABLE 7.1. Modeled Stratospheric Aerosol Extinction 



MOD 4 

MOD 5 

MOD 6 

X (tun) 

EXT, m" 1 

EXT, m" 1 

EXT, rrf 1 

.385 

.314 x 10~® 

1.95 x 10~® 

8.86 x 10“ 6 

.448 

.292 x 10~® 

1.66 x 10"® 

8.69 x 10 -6 

.453 

.291 x 10 _ ® 

1.65 x 10~® 

8.71 x 10 -6 

.525 

.200 x 10*® 

1.43 x 10" 6 

9.06 x 10 -6 

.600 

.226 x 10“® 

1.248 x 10" 6 

9.20 x 10* 6 

.940 

.107 x 10~® 

.715 x 10" 6 

7.82 x 10 -6 

1.02 

.089 x 10' s 

.631 x 10~ 6 

7.40 x 10" 6 

m (kg m -9 ) 

6.5 x 10“ U 

5.01 x lO" 10 

4.09 x 10" 9 

N (m“ 3 ) 

.117 x 10 7 

1 

.273 x 10® 

.155 x 10® 
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The results of this analysis are shown in Table 7.1, together with the 
aerosol mass density (N) and particle concentation (n). Points to note are: 

1. The large effect on the extinction caused by the volcanic eruptions, 
particularly at the longer wavelengths; 

2. The strong wavelength variation in extinction for the prevolcanic 
model and the lack of variation in the post El Chichon model. 

Both these results reflect the effects of the larger aerosols present 
after the eruptions, particularly after that of El Chichon. 

7.3 High-Altitude Lidar Modeling 

Recent flights of the NASA Shuttle spacecraft have shown that the 
re-entry region between 70 and 100 km altitude is one in which large 
molecular density fluctuations occur (e.g., Price 1983). These fluctuations, 
which may be as great as -30% to +60% of the ambient density, are probably 
associated with the passage through this region of atmospheric gravity 
waves. The potential value of being able to study these fluctuations from 
the ground has led to an investigation of the possibility of upgrading and 
using the NASA-LaRC 48" lidar system for this purpose. Modeling studies of 
the theoretical performance of such an upgraded system have been carried 
out. Some of the results of this modeling are given in Tables 7.2, 7.3, and 

7.4 which show the expected percentage uncertainty in the molecular 
backseat tered signal from different altitudes for a 10-minute observational 
period. In all the cases shown here, it has been assumed that photon 
counting techniques are used and that the noise depends only upon signal 
fluctuations and background lighting conditions. Reasonable figures are used 
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TABLE 7.2. 


Theoretical Performance Figures for 48" Lidar 
System — Nighttime, Interference Filter 


(A) SYSTEM PARAMETERS 

2 

Receiving Mirror Area - 1.00 m 

Full Angle Beam Divergence - 0.200E-03 r 

Wavelength - 0.5300E-06 m 

Overall Efficiency - 0.08 

Pulse Energy - 0.4 J 

Firing Rate - 10 Hz 

Filter Bandwidth - 0.100E-08m 

(B) ATMOSPHERIC PARAMETERS 

Atmospheric Transmission - 0.70 
Sky Brightness - 0.70E+00 R.A.”^ 

(C) OPERATIONAL PARAMETERS 

Height Increment - 2.0 km 
Integration Time - 600 S 

System Constant - 8.36E+19 


Height 

(km) 


Signal 

Count 


Noise % Uncertainty 

Count (10 Minutes Integration) 


65 

4.26 

1.85 

X 

10 4 

0.63% 

70 

1.90 

1.85 

X 

io-4 

0.94% 

75 

0.815 

1.85 

X 

io-4 

1.43% 

80 

0.331 

1.85 

X 

io-4 

2.24% 

85 

0.121 

1.85 

X 

io-4 

3.72% 

90 

0.0428 

1.85 

X 

10-4 

6.25% 

95 

0.0150 

1.85 

X 

10-4 

10.59% 

100 

0.0057 

1.85 

X 

10-4 

17.40% 

105 

0.0022 

1.85 

X 

10-4 

28.56% 

110 

0.0009 

1.85 

X 

10-4 

45.93% 


10 min - 10% Uncertainty Altitude = 94 km 



TABLE 7.3. Theoretical Performance Figures for 

48" Lidar System — Daytime, Interference 
Filter 


(A) SYSTEM PARAMETERS 

2 

Receiving Mirror Area - 1.00 m 

Full Angle Beam Divergence - 0.200E-03 r 

Wavelength - 0.5300E-06 m 

Overall Efficiency - 0.08 

Pulse Energy - 0.4 J 

Firing Rate - 10 Hz 

Filter Bandwidth - 0.100E-08 m 

(B) ATMOSPHERIC PARAMETERS 

Atmospheric Transmission - 0.70 
Sky Brightness - 8.0E+6 

(C) OPERATIONAL PARAMETERS 

Height Increment - 2.0 km 
Integration Time - 600 S 

System Constant - 8.36 E+19 


Height 

(km) 

Signal 

Count 

Noise 

Count 

% Uncertainty 
(10 Minutes Integration) 

65 

4.26 

2119 

13.9% 

70 

1.90 

2119 

31.2% 

75 

0.82 

2119 

72.9% 

80 

0.33 

2119 

179.6% 

85 

90 

95 

100 

105 

110 

10 min - 

10% Uncertainty Altitude 

= 63 km 
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TABLE 7.4. Theoretical Performance Figures for 

48" Lidar System — Daytime, Interference 
Filter Plus Fabry-Perot Etalon 


(A) SYSTEM PARAMETERS 

2 

Receiving Mirror Area - 1.00 m 

Full Angle Beam Divergence - 0.200E-03 r 

Wavelength - 0.5300E-06 m 

Overall Efficiency - 0.04 

Pulse Energy - 0.4 J 

Firing Rate - 10 Hz 

Filter Bandwidth - 1.0E-11 m (0.1 A) 

(B) ATMOSPHERIC PARAMETERS 

Atmospheric Transmission - 0.70 
Sky Brightness - 8.0E+6 R.A“1 

(C) OPERATIONAL PARAMETERS 

Height Increment - 2.0 km 
Integration Time - 600 S 

System Constant - 4.18E+19 


Height 

Signal 

Noise 

% Uncertainty 

(km) 

Count 

Count 

(10 Minutes Integration) 

65 

2.13 

10.66 

2.17% 

70 

0.95 

10.66 

4.83% 

75 

.0.407 

10.66 

10.55% 

80 

0.166 

10.66 

25.67% 

85 

0.0603 

10.66 

70.13% 


90 

95 

100 

105 

110 


10 min - 10% Uncertainty Altitude - 75 km 
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for angular beam-widths, system efficiencies and sky brightness; no allowance 
has been made for any post-detector noise or signal loss. 

The conditions assumed in the tables are: 

Table 7.2. Night-time with 10 A interference filter 

Table 7.3. Day-time with 10 A interference filter 

Table 7.4. Day-time with Fabry-Perot Etalon 

In Table 7.2, the signal-to-noise ratio is determined mainly by signal 
fluctuations and useful measurements are obtainable up to about 90 km. In 
Table 7.3, for day-time, the sky brightness is approximately 10^ greater than 
at night and this severely limits the system performance for altitudes above 
65 km. Table 7.4 shows the improvement that may be expected if a 0. 1 A 
Fabry-Perot Etalon is used in conjunction with the interference filter. 
Considerable improvement in system performance occurs; it now being possible 
to obtain useful information up to about 75 km. 

Results of further calculations, including modeling frequency doubled 
Nd-YAG, and ruby laser transmitters with different laser output energy levels 
have been presented in the final report to NASA, Contract NAS 1-17959 (Kent 
and Wang, 1986). The latter report has also included an examination of the 
science requirements for possible density fluctuation studies. 
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8. CONCLUSIONS 


A detailed study has been made of the SAM II data set from October 1978 
to February 1983 and of the entire SAGE I data set from February 1979 to 
November 1981. The emphasis of this study has been on the validation of these 
data sets by their use in various analyses of the dynamical, physical and 
chemical processes taking place in the stratosphere. Much of the analysis 
has been climatological; this has included studies of the statistics of polar 
stratospheric clouds and of the global scale changes taking place in the 
stratosphere following volcanic injection. Related analyses have included 
the preparation of plots, tables and films showing these climatologies in as 
direct and simple a manner as possible. In the course of this work, various 
anomalies and inconsistencies in the data have been noted, this information 
has either been given to NASA to permit updating of successive editions of 
the archived data or has been specifically listed for user application (e.g.. 
Table 2-6). 

The SAGE I and SAM II data have also been used for the study of 
atmospheric phenomena, such as the changes in aerosol and ozone 
concentrations taking place in the vicinity of the north polar vortex. The 
SAGE I multiwavelength data have been used to investigate changes in the 
aerosol size distribution and the ozone data to study winter-time eddy fluxes 
in relation to planetary wave activity. As part of the direct data 
validation program, SAGE I and SAM II, 1.0 um aerosol data have been 
intercompared wherever the satellite events have occurred close together in 
space and time. 

These studies have shown the quality of the data to be excellent. 
Results of the scientific analyses have been presented at scientific 
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conferences and have been published in relevant journals. The data, now 
archived, will represent an invaluable information source for many years, the 
analyses described in this report representing only a fraction of its 
possible uses. 
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APPENDIX A 


A COMPARATIVE STUDY OF AEROSOL EXTINCTION MEASUREMENTS MADE 
BY THE SAM II AND SAGE SATELLITE EXPERIMENTS 


by 


Glenn K. Yue 1 

Institute for Atmospheric Optics and Remote Sensing (IFAORS) 

M. P. McCormick and W. P. Chu 
NASA Langley Research Center 


Abstract 


SAM II and SAGE are two satellite experiments designed to measure 
stratospheric aerosol extinction using the technique of solar occultation or 
limb extinction. Although each sensor is mounted aboard a different 
satellite, there are occasions when their measurement locations are nearly 
coincident, thereby providing opportunities for a measurement comparison. In 
this paper, the aerosol extinction profiles and daily contour plots for some 
of these events in 1979 are reported. The comparisons shown in this paper 
demonstrate that SAM II and SAGE are producing similar aerosol extinction 
profiles within their measurement errors and that since SAM II has been 
previously validated, these results show the validity of the SAGE aerosol 
measurements. 


(From J. of G eo phy. Res ., 89 , 1984, 5321-5327.) 
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APPENDIX B 


POLAR STRATOSPHERIC CLOUD SIGHTINGS BY SAM II: 1978-1982 


by 


U. 0. Farrukh 

Institute for Atmospheric Optics and Remote Sensing, 
Hampton, Virginia 23666 


Abstract 


Satellite data sets, obtained by the SAM II experiment, are analyzed to 
study the characteristics and distribution of Polar Stratospheric Clouds 
(PSC). Mappings and tables were made of their distribution and probability of 
formation. A 50% probability of formation was found to exist around a minimum 
stratospheric temperature of 193°K. The highest probability of cloud 
formation was found to be close to 0° longitude. However, some differences in 
clouds altitude distribution were found to exist between the northern and 
southern polar regions. Satellite observations inside and outside the 
boundaries of the asymmetric northern polar night jet are made daily. It was 
possible to conclude that these clouds are formed only inside the polar jet 
boundaries in the northern hemisphere. In general, the number of PSC events 
that were observed in the southern hemisphere were three to four times the 
ones observed in the northern one. 


(From NASA Contractor Report 177995, March 1986.) 
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APPENDIX C 


CHARACTERISTICS OF POLAR STRATOSPHERIC CLOUDS 
AS OBSERVED BY SAM II, SAGE, AND LIDAR 


by 


M. P. McCormick 

NASA Langley Research Center Hampton, Virginia 
Patrick Hamill 

Systems and Applied Sciences Corporation Palo Alto, California 

and 

U. 0. Farrukh 

Institute for Atmospheric Optics and Remote Sensing (IFAORS) 

Hampton, Virginia 


Abstract 


The discovery of polar stratospheric clouds (PSC's) in the Arctic and 
Antarctic stratosphere during winter is described, and the directly observed 
and implied properties of these clouds are discussed. It is proposed that the 
more familiar "mother-of-pearl" or "nacreous" clouds are a special subset of 
PSC's. The size, location, prevalence and temperature dependence of the 
clouds as measured by the SAM II and SAGE satellite systems are outlined. 
Airborne lidar measurements have recently demonstrated that the PSC phenomenon 
is most probably associated with an extended stratospheric cloud bank existing 
within the cold polar vortex region during the winter period with the PSC's 
bounded by a 188 temperature isotherm. The PSC's probably exist at a 50 
percent frequency within the 193 K isotherm. Using the observed information 
on the cloud extinction and change in location with time we consider possible 
formation mechanisms, the size of the cloud particles, and show the descending 
motion of the cloud during wintertime. 


(Reproduced from the Journal of the Meteorological Society of Japan , 1985, 
63, 267-276.) 
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APPENDIX D 


SPATIAL CHANGES IN THE STRATOSPHERIC AEROSOL ASSOCIATED WITH THE 

NORTH POLAR VORTEX 


by 


M. P. McCormick 
NASA Langley Research Center 

C. R. Trepte 
SASC Technologies, Inc., 

and 

G. S. Kent 

Institute for Atmospheric Optics and Remote Sensing (IFAORS) 

Hampton, Virginia 


Abstract 


In late January and early February 1983, observations made by the 
Stratospheric Aerosol Measurement (SAM II) satellite system showed that 
aerosol extinction profiles measured within the northern polar vortex differed 
significantly above 18 km from those measured outside the vortex. Values of 
the calculated optical depths above 18 km for February 1, 1983, are lower by 
approximately one order of magnitude within the polar vortex than those 
outside. Similar differences were found in the aerosol back-scattering 
profiles obtained using an airborne lidar system when crossing the polar 
vortex. Since potential vorticity at a constant altitude is not conserved 
across the polar vortex, horizontal adiabatic transport does not occur. 


(Reproduced from Geo ph ys. Res. Lett. , 1983, 941-944.) 
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APPENDIX E 


VARIATION IN THE STRATOSPHERIC AEROSOL ASSOCIATED WITH 
THE NORTH CYCLONIC POLAR VORTEX AS MEASURED 
BY THE SAM II SATELLITE SENSOR 


by 


G. S. Kent 

Institute for Atmospheric Optics and Remote Sensing (IFAORS) 

Hampton, Virginia 


C. R. Trepte 

SASC Technologies, Inc., Hampton, Virginia 23666 


U. 0. Farrukh 

Institute for Atmospheric Optics and Remote Sensing (IFAORS) 
Hampton, Virginia 


and 


M. P. McCormick 

NASA Langley Research Center, Hampton, Virginia 


Abstract 


Aerosol extinction data obtained by the Stratospheric Aerosol Measurement 
II (SAM II) satellite instrument during the 1979/1980 Northern Hemisphere 
winter season have been analyzed in relation to the cyclonic polar vortex. A 
Synoptic approach has been employed to study the behavior of aerosol 
extinction ratio and optical depth between altitudes of 8 and 30 k as a 
tracer of mean atmospheric motions in and near the polar vortex. As the polar 
vortex intensifies, a gradient of extinction ratio is established across the 
polar-night jet stream, which is associated with subsidence within the vortex. 
Maximum subsidence occurs at the center of the vortex. Calculated descent 
rates relative to isentropic surfaces are of the order of 8 x 10“Sn s _i near 
20 km, at the center of the vortex between September and December. Below an 
altitude of 14 km, taken as the base of the vortex, and outside the vortex, 
horizontal movements occur freely, masking any systematic vertical motions. 
Extinction enhancements by polar stratospheric clouds and changes produced by 
sudden warmings in the second half of winter have prevented a similar study 
for this period. An estimate of the aerosol mass transferred downward through 
the base of the vortex for the entire season is 7000 tonnes. Comparison of 
the inferred stratospheric motions with earlier studies using radioactive 
tracers show good agreement. 


(Reproduced from J. Atmos. Sc i., 42^ 1985, 1536-1551.) 
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APPENDIX F 


LATITUDINAL AND ALTITUDINAL VARIATION OF SIZE DISTRIBUTION 
OF STRATOSPHERIC AEROSOLS INFERRED FROM SAGE AEROSOL 
EXTINCTION COEFFICIENT MEASUREMENTS AT TWO WAVELENGTHS 


by 


Glenn K. Yue 1 and Adarsh Deepak 
Institute for Atmospheric Optics and Remote Sensing (IFAORS) 

Hampton, Virginia 

Abstract 

A method of retrieving aerosol size distribution from the measured 
extinction of solar radiation at wavelengths 0.45 pm and 1.0 pm has recently 
been proposed. This method is utilized to obtain latitudinal and altitudinal 
variations of size distributions of stratospheric aerosols from the SAGE data 
for March 1979. Small particles are found in the lower stratosphere of the 
tropical region, and large particles are found at higher altitudes and 
latitudes in both hemispheres. Results of this study are consistent with the 
suggestion that the upper troposphere in tropical regions is a source of 
condensation nuclei in the stratosphere, and they become mature as they move 
to higher altitudes and latitudes. 


(From Geoph ys . Res ♦ Lett . , 1 1 , 1984, 999-1002.) 
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APPENDIX G 


SAGE AND SAM II MEASUREMENTS OF GLOBAL 
STRATOSPHERIC AEROSOL OPTICAL DEPTH 
AND MASS LOADING 


by 


G. S. Kent 

Science and Technology Corporation 
Hampton, Virginia 

and 

M. P. McCormick 
NASA Langley Research Center 
Hampton, Virginia 


Abstract 


Several volcanic eruptions between November 1979 and April 1981 have 
injected material into the stratosphere. The SAGE and SAM II satellite 
systems have measured, with global coverage, the 1-pm extinction produced by 
this material and examples of the data product are shown in the form of global 
maps of stratospheric optical depth and altitude-latitude plots of zonal mean 
extinction. These data, and that for the volcanically quiet period in early 
1979, have been used to determine the changes in the total stratospheric mass 
loading. Estimates have also been made of the contribution to the total 
aerosol mass from each eruption. It has been found that between early 1979 
and mid- 198 1, the total stratospheric aerosol mass increased from a background 
level of approximately 570,000 metric tons to a peak of approximately 
1,300,000 metric tons. 


(From J. Geophys. Res ., 89 , 1984, 5303-5314.) 
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APPENDIX H 


BEHAVIOR OF ZONAL MEAN AEROSOL EXTINCTION RATIO AND ITS 
RELATIONSHIP WITH ZONAL MEAN TEMPERATURE DURING 
THE WINTER 1978-1979 STRATOSPHERIC WARMING 


by 


Pi-Huan Wang 

Institute for Atmospheric Optics and Remote Sensing (IFAORS) 

Hampton, Virginia 

and 

M. P. McCormick 
NASA Langley Research Center 
Hampton, Virginia 


Abstract 


The behavior of the zonal mean extinction ratio in the lower stratosphere 
near 75 N and its relationship with the zonal mean temperature during the 
January-February 1979 stratospheric sudden warming have been investigated 
based on the satellite sensor SAM II (Stratospheric Aerosol Measurement) and 
auxiliary meteorological measurements. The results indicate that distinct 
changes in the zonal mean aerosol extinction ratio occurred during this 
stratospheric sudden warming. It is also found that horizontal eddy transport 
due to planetary waves may have played a significant role in determining the 
distribution of the zonal mean aerosol extinction ratio. 


(From J. of Geophys. Res. , 90, 1985, 2360-2364.) 
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APPENDIX I 


VARIATIONS IN STRATOSPHERIC AEROSOL OPTICAL DEPTH DURING 

NORTHERN WARMINGS 


by 


Pi-Huan Wang 

Institute for Atmospheric Optics and Remote Sensing (IFAORS) 
Hampton, Virginia 

and 

M. P. McCormick 
NASA Langley Research Center 
Hampton, Virginia 


Abstract 


In this paper the properties of the stratospheric aerosol optical depth 
(above 50 mbar) have been studied by using aerosol extinction profiles (at 1 
pm) derived from the Stratospheric Aerosol Measurement (SAM II) and 
Stratospheric Aerosol and Gas Experiment (SAGE) during warming periods in the 
northern hemisphere. It is shown that, during the disturbed periods in 
winter, low values of aerosol optical depth (<0.2 x 10-3) are found within the 
low-pressure system(s) (at the 30-mbar pressure surface), while high values 
are found outside. Similar characteristics are found to exist for the 
simultaneously observed SAGE O 3 and NO 2 columnar density distributions. 
Strong longitudinal gradients are shown with the low values within and 
wherever the vortex exists. This characteristics is maintained during and 
after the circumpolar vortex is disturbed, even after breakdown, indicating an 
isolation of the material within the vortex. 


(From J. Geophys. Res ., 90 , 1985, 10,597-10,606.) 
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APPENDIX J 


retrieval of stratsospheric aerosol size distribution from 

ATMOSPHERIC EXTINCTION OF SOLAR RADIATION AT 
TWO WAVELENGTHS 

by 

Glenn K. Yue and Adarsh Deepak 
Institute for Atmospheric Optics and Remote Sensing (IFAORS) 


Abstract 


The possibility of retrieving size distribution of aerosol particles in 
the stratosphere by measuring the extinction of solar radiation traversing the 
aerosol medium at two different wavelengths is explored in this work. This 
paper presents a parametric study of the effects of size distribution and 
composition of stratospheric aerosols on the value R, the ratio of extinction 
of solar radiation at wavelength 0.45 pm to the extinction at wavelength 1.0 
urn. The aerosol size distributions under study are nine analytical 
expressions including most stratospheric aerosol models used by investigators. 
The aerosol compositions under study are the supercooled sulfuric-acid 
droplets, with different weight percentages of ^SO ^ in the aerosol. It is 
found that R is not very sensitive to either the composition or the radii 
limits of stratospheric aerosols under consideration but is quite sensitive to 
the value of the variable parameter governing the mode radius of the size 
distribution of aerosol particles. Based on the results of this parametric 
study, a method of retrieving the size distribution of aerosol particles in 
the stratosphere from the experimental results of the extinction of solar 
radiation at two wavelengths is proposed. 


(From Appl . Opt ., £2, 1983, 1639-1645.) 
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